Transport Properties of Multiple Quantum Dots Arranged in Parallel: Results from 

the Bethe Ansatz 
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In this paper we analyze transport through a double dot system connected to two external leads. 
Imagining each dot possessing a single active level, we model the system through a generalization of 
the Anderson model. We argue that this model is exactly solvable when certain constraints are placed 
upon the dot Coulomb charging energy, the dot-lead hybridization, and the value of the applied gate 
voltage. Using this exact solvability, we access the zero temperature linear response conductance 
both in and out of the presence of a Zeeman field. We are also able to study the finite temperature 
linear response conductance. We focus on universal behaviour and identify three primary features 
in the transport of the dots: i) a so-called RKKY Kondo effect; ii) a standard Kondo effect; and 
iii) interference phenomena leading to sharp variations in the conductance including conductance 
zeros. We are able to use the exact solvability of the dot model to characterize these phenomena 
quantitatively. While here we primarily consider a double dot system, the approach adopted applies 
equally well to N-dot systems. 



PACS numbers: 72.15.Qm, 73.23.Hk, 85.35.Gv, 71.10.Pm 



I. INTRODUCTION 

Quantum dot devices exhibit a wide range of strongly 
correlated phenomena. In the simplest of cases, single dot 
devices exhibit classical Kondo physics the archetype 
of strong correlations in impurity systems. This physics 
arises in dots due to their natural localization of elec- 
tronic (and hence spin) degrees of freedom. Measuring 
transport properties provides the most effective means to 
probe the Kondo physics in such dots. The formation of 
the Abrikosov-Suhl resonance, a dynamical enhancement 
in the density of states at the Fermi energy, facilitates 
transport through the dot from one lead to the other. 
Manipulations of this resonance through adjustments of 
temperature, bias, or external fields are readily observed 
through monitoring the device's conductance. 

At the next level of complexity are single level dots 
where multiple energy levels are relevant to transport. 
Much attention has been paid to the situation where two 
levels in a single dot are anomalously close In such 

a case the dot can be tuned through a singlet-triplet tran- 
sition. On the triplet side a two stage Kondo effect can be 
realizcd££ Here two channels of electrons coupled with 
different strengths to the dot triplet successively reduce 
the impurity from S = 1 to S 1 = 1/2 to <S = in a two 
stage process. 

Quantum dot devices with multiple active dots offer 
the same range of phenomenology as single multi-level 
dots but with greater versatility as level spacing and in- 
teractions are not precisely controllable within a single 
dot. In particular multiple quantum dot devices offer the 
opportunity to study non-trivial mixtures of strongly cor- 
related physics and interference phenomena. Any Kondo 
physics that arises in such systems will be strongly af- 
fected by the presence of multiple tunneling paths af- 
forded by the multiple dots. Even in the absence of in- 
teractions, multiple tunneling paths lead to richly struc- 



tured transport behaviour predicated upon interference 
effects. 

Multiple dot devices come in a number of guises. In 
one instance they are manufactured from gated semi- 
conducting heterostructuresi^iSiiiS An important feature 
that the semiconducting multi-dots share with their sin- 
gle dot cousins is their tunability. By adjusting voltages 
applied to various gates which define the dots, it is possi- 
ble to change fundamental parameters of the dot system 
allowing experimentally exploration of all of the dots' 
various physical regimes. In a second instance they can 
be constructed from carbon nanotubesi 11 ! 12 With appro- 
priate gating, carbon nanotube dots can also be made 
fully tunable i 11 ' 12 Much of the experimental interest in 
multiple quantum dots comes from their possible realiza- 
tion as elementary qubits in a quantum computer j 10 ' 13 ' 14 
As the spin degrees of freedom on the dot serve as the 
states of the qubit, the question of spin decoherence due 
to coupling with the nuclear background becomes an im- 



portant question 
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Carbon nanotube dots have 



the advantage that their spin decoherence times are ex- 
pected to be longer due to a weaker nuclear background 
in carbon than in GaAsi n i 12 

The promised rich phenomenology and the ability to 
realize multiple dot devices together have led to a flurry 
of theoretical work in this area. Combinations of inter- 
ference effects and Kondo physics in parallel dots have 
been studied in many works i 16 i 17 ' 18 i 19 ' 20 i 21 ' 22 In response 
to the observations of a competition between an effec- 
tive Ruderman-Kittcl-Kasuya-Yosida (RKKY) coupling 
and a Kondo effect in a pair of double dots in series^ a 
number of works have studied the physics of RKKY ef- 
fects in quantum dotS i 17 ' 23 i 25 Techniques that have been 
used to study multiple dots include infinite U slave boson 
mean field theory , 16 ' 17 i 18 ' 19 i 24 ' 25 i 26 i 27 numerical renormal- 
ization group ) 16 ' 21 i 24 ' 28 i 29 exact diagonalization ) 20 i 30 ap- 
peal to asymptotic limits^ i 5 ' 23 i 31 , a perturbative renor- 
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malization group 3 -^, and a functional renormalization 
group 33 ' 34 . 

In this work we present a distinct approach to ana- 
lyzing such systems: we study such systems under the 
rubric of integrability. Like the numerical renormaliza- 
tion group it offers exact results but with the advantage 
of analytic control. Integrability has long been known 
to be able to access thermodynamic properties of single 
level dot systems^ More recently it has been shown to 
be successful in computing equilibrium transport proper- 
ties of single level quantum dots^ In particular, these 
techniques were able to compute the finite temperature 
linear response conductance over a number of decades in 
temperature in good agreement with the numerical renor- 
malization group. The aim of this paper is to apply this 
approach to multiple quantum dots. It is an extension 
and elaboration of work first reported in Ref. (f37h . 

To this end, we examine a generalized Anderson model 
coupling two leads (I = 1, 2) to an array of dots in paral- 
lel: 

/OO 
dxc\ a d x ci„ + ^2 Vi a (c\ aa d aa + h.c) 

la -°° era 

+ / ] ^dg^gg + 7J Ugg'Tltania. (1.1) 
era aa' 

Here the ci a /d aa specify electrons living in the leads and 
the dots, a indexes the various levels on the dots in the 
system. Vi a measures the tunneling strength between the 
dot level a and lead I. U aa > characterizes the Coulom- 
bic repulsion between electrons of opposite spin living in 
levels a and a'. A schematic of this Hamiltonian for two 
dots, i.e. a = 1, 2, in given in Figure 1. 



Ugc 



U\ F a =T a r J Eda ~ £da' — £d- (1-3) 




FIG. 1: A schematic of two dots in parallel. 

We will argue in this paper that the above Anderson- 
like Hamiltonian is exactly solvable in two cases: 

Vi a / Vac = Vic/Vw, 

U a + 2e da = U a . +2e da ,; (1.2) 



and 



V la /V 2a = V la ,/V 2 g, 



The first condition in both cases guarantees that 
only a single effective channel in the leads couples 
to the dots. This condition is commonplace in the 
literature^L^^^^^i 3 ^^^. The second condition 
implies that for fixed Coulomb repulsion of the dots, the 
bare level separation of the two dots, |e<ji — e^l, is also 
fixed. This still permits the dot levels to be moved in 
unison, say by the application of a uniform gate voltage. 
The final condition amounts to the insistence that the en- 
ergy scale associated with the charge degrees of freedom, 
i.e. sfUjTi, must be the same on both dots. However 
this is not a serious constraint as we are much more in- 
terested in physics occurring on a scale associated with 
spin fluctuations which are not so constrained. 

In the bulk of the paper, we focus this approach on 
studying double dots (i.e. a = 1,2 in the above model). 
Wc identify a number of interesting phenomena. At the 
particle-hole symmetric point of the double dot system 
we find that the ground state of the system is a singlet 
and is characterized by a Kondo-like resonance. This im- 
plies that the RKKY interaction between the dots acts in 
an unusual fashion. In the standard picture, one expects 
cither that an antifcrromagnctic RKKY interaction binds 
the two electrons into a singlet but with any Kondo effect 
absent or that a ferromagnetic RKKY interaction binds 
the electrons into a triplet which is subsequently reduced 
upon Kondo screening to a doublet. But here we find nei- 
ther scenario. Instead wc find an effective antifcrromag- 
netic RKKY interaction that cooperates, not competes, 
at an equivalent energy scale with Kondo singlet forma- 
tion. Wc term this an RKKY Kondo effect. This RKKY 
effect is clearly much different than the pcrturbative fer- 
romagnetic interaction^ between dots appropriate for a 
high temperature analysis of the dot system.— Using the 
Bethe ansatz, we are able to describe quantitatively the 
associated Abrikosov-Suhl resonance together with how 
it evolves under applications of finite gate voltage (thus 
destroying particle-hole symmetry), finite magnetic field, 
and finite temperature. 

We also identify a phenomena that we believe is under- 
appreciated. At values of the gate voltage where the dot 
system is transitioning from the mixed valence regime to 
the empty orbital regime, the number of electrons dis- 
placed, ndis, in connecting the leads to the dots is neg- 
ative, ridis is a sum of both the localized electrons on 
the dots and changes in the electron density in the leads. 
Our analysis allows us to argue that this second contri- 
bution to ridis is finite (and negative). This phenomena 
is associated in general with the presence of interference 
and for example has also been predicted to occur in dots 
exhibiting Fano resonances^ 

A detailed outline of the paper is as follows. It is di- 
vided in two overall parts. The first part, comprising 
Sections II through IV, is technical in nature. In Sec- 
tion II we both introduce the model of N-dots in parallel 
and demonstrate that it admits an exact solution via the 
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Bethe ansatz. To do so we construct both the one and 
two particle wavefunctions explicitly. This allows the ex- 
traction of the S-matrix and thus the ability to construct 
implicitly the N-particle wavefunctions through the al- 
gebraic Bethe ansatz. We also sketch briefly the explicit 
construction of the N-particle wavefunctions. This is nec- 
essary to ensure the posited Bethe ansatz wavefunctions 
are indeed eigenfunctions of the dot Hamiltonian. Finally 
in Section II we provide a general analysis of the Bethe 
ansatz equations in the continuum limit both at zero and 
finite temperature. 

In Section III wc sketch how the transport properties 
of the double dot can be accessed via the Bethe ansatz 
equations. We do so on the same basis that the mag- 
netoconductance was computed in the Kondo model. 42 
There the scattering phase of an electron off the spin im- 
purity was computed by analyzing the 1/L (where L is 
the system size) corrections of the Bethe ansatz equa- 
tions. Here the analysis is more involved as Anderson- 
type models have non-trivial excitations involving both 
charge and spin degrees of freedom (unlike the Kondo 
model where the charge sector decouples). We, however, 
are still able to relate the scattering phase to an impu- 
rity density of states. As this latter quantity is readily 
extracted from the Bethe ansatz equations, we are in po- 
sition to characterize the transport through the dot. In 
Section III we sketch how this proceeds both at zero and 
finite temperature. 

Up to this point in the article, we have made no at- 
tempt to provide specific solutions to the integral equa- 
tions arising out of the analysis of the Bethe ansatz 
quantization condition. In Section IV of the paper this 
changes. We both specialize to the case of two dots 
and wc outline how the equations describing transport 
through two-dot systems can be solved analytically in 
the two particular limits of large and small bare level 
separation, i.e. \e dl - e d %\ > V? 2 and \e dl - e d2 | < V? 2 . 
This ends the technical portion of the paper. 

In the second part of the paper, Sections V through 
VII, wc discuss the various features of transport arising 
in a double parallel dot system. These sections can be 
read (mostly) independently from the Bethe ansatz anal- 
ysis in the paper's first part. In Section V wc study zero 
temperature transport both in and out of an applied Zee- 
man field. In zero field we study the conductance through 
the double dots as a function of gate voltage. We iden- 
tify three regions of gate voltage: i) near the particle-hole 
symmetric point where two electrons sit on the dots and 
where what we term an RKKY Kondo effect is present; 
ii) a region in gate voltage where one electron sits on the 
dots and a standard Kondo effect takes place; and iii) 
a region in gate voltage marking the transition between 
a mixed valence and empty orbital regime where marked 
interference effects are present. We then study how these 
various phenomena are affected by the introduction of a 
Zeeman field. In Section VI we study the finite tem- 
perature linear response conductance at the particle-hole 
symmetric point (i.e. how the introduction of a tempera- 



ture destroys the RKKY Kondo effect). In Section VII we 
end the paper considering two items. We first examine 
how the Friedel sum rule operates in the region where 
strong interference effects are present. And second we 
consider how breaking the integrability constraints, i.e. 
Eqns. (|1.2|) and (|f .3|) . affects the physics we have dis- 
cussed in the manuscript. On this last point we conclude 
that the effects are in general perturbativc, i.e. small vi- 
olations of the constraints leads only to small changes in 
the physics. 



II. EXACT SOLVABILITY OF N-DOTS IN 
PARALLEL 

A. Integrability of Model: Bethe Ansatz Equations 

To begin to analyze the model in Eqn. (1.1) we first 
map the problem to an Anderson model involving a single 
effective lead. To do so we need to assume the ratio 
of left/right lead couplings are equal, i.e. V\ a /V2 a = 
Vi a >/V 2a >. Writing c e/o = {V 1/2a ci ± V 2/la c 2 )/\/2T a , 
with r a = (Vf a + V 2a )/2, the Hamiltonian factorizes into 
an even and an odd sector: 

/oo 
dx c\ a d x c e(7 + ^2 V / 2T^(4 CTQ d CJQ + h.c) 

ia 00 era 
era aa' 

/oo 
dx c\ a d x c oa . 
OO 



(2.1) 



la 



Only the even Hamiltonian, coupling directly to the dot 
degrees of freedom, is non-trivial. 

To determine under what conditions 7i e is exactly 
solvable by Bethe ansatz we compute both the one and 
two electron eigenstates. The one particle wave function 
takes the form, 



dx{g a (x)cl(x)} + e acr d) aa . 



|0). (2.2) 



We solve the Schrodinger equation, H. e \^p) = an d 
find g(x) is of the form, 



g a (x) = 9{x)e lqx+l ^ L 



:«(<!) 



(2.3) 

with S(q) the scattering phase of the electron off the dot 
array equaling, 



%) = -2tan" 1 (^ 



(2.4) 



The total scattering phase is a sum of the arguments of 
tan -1 of scattering phases off the individual dots in the 
system. 

The two particle eigenfunction in the S z = sector 
takes the form 



/oo poo 
dxidx 2 g(xi 1 X2)c^(x 1 )c^(x 2 ) + / dx 
-oo „, J — oo 



[e a {x){c]{x)d} ai - c\(x)d\J]+Y^f aa 'd\ a dl a , 

OLOt' 



|0). (2.5) 



Again solving the Schrodinger equation Ti e \ip) = (q + 
p)\i/j) gives g(xi,x 2 ) to be 

g(x 1 ,x 2 ) = g q {xi)g p (x 2 )(j)(xi - x 2 ) 

+g P (xi)g q (x 2 )4>(x 2 - xi). (2.6) 

Here g q / p (x) and e q / p are the coefficients appearing in 
the one particle wavefunction corresponding to energies 
q/p. </>(x) governs the scattering when two electrons are 
interchanged. It takes the form <p{x) = l+i"/(q,p)sign(x). 
We find that j(q,p) must satisfy for arbitrary a, 



j{q,p) 



1 1 



E 



p a q 



Pa>q 



2q-p^-f p a ,q D 

a' 

2T a / U aa > 



q-p 



(2.7) 



where p a /q a = p - e da /q - e da . For consistency, j(p, q) 
must be independent of a particular vale of a. This only 
happens if one of two sets of conditions holds: 



case i: 



U a + 2e da = U a ' + 2e da r : 



case n: 



Uaa' 



(~<Io 



(2.8) 



(2.9) 



Case i) would be appropriate to describing well separated 
single level dots while case ii) would be appropriate to 
describing a single dot with degenerate levels. 

Exact solvability is predicated on how "f(q,p) deter- 
mines the scattering matrix of the two electrons. The 
scattering matrix has the general spin (SU(2)) invariant 
form 



Sl\ll = b{p, q)fa\%l + c(p, q)Palai, 



(2.10) 



where here {<Ji,a 2 }/{aii,a 2 i} represent the incoming 

and outgoing spins. Here Ialal = $cr\5al is the identity 

matrix while Palal = ^a\^a\ is the permutation matrix. 
The coefficients, b(p,q) and c(p,q) 7 are determined by 
j(p, q) from the two relations 



b(p, q) - c(p, q) 
b(p, q) + c(p, q) 



(xi - x 2 ) 



1 + 47 (p, q) 



(j){x 2 -x 1 ) l-ij(p,q) 

1. (2.11) 



The first relation results from the two particle eigenfunc- 
tion in the S- = sector solved above while the second 



time 





'2 °3 



FIG. 2: A graphical representation of the Yang-Baxter equa- 
tion. Three particles cti,o"2, and 03 scatter into a'\,a'i, and 
cr" in two different ways. Under the Yang-Baxter relation, 
the two scattering processes are equivalent. 



relation arises from the corresponding eigenfunction in 
the S z = ±1 sector where interactions between spins are 
trivial. 

In order for the Hamiltonian Ti e to be integrable, the 
above S-matrix must satisfy the Yang-Baxter relation. 
The Yang-Baxter relation governs the scattering of three 
electrons, u\,(J 2l and 03. It enforces the equivalency of 
different scattering orders. a\ may first scatter off a 2 and 
then <7 3 followed by a scattering of a 2 and 03 or a 2 and 
a 3 may scatter first followed by successive scattering of 
<j\ off 02 and 03. This is graphically encoded in Figure 
2. Algebraically it is given by 

Sl\ll (pi,p 2 ) S a jll (pi,p 3 ) (p 2 ,p 3 ) = 



l(P2,P3)SZl (Pi.Ps)^ (Pi,pa).(2.12) 



If this relation holds, the scattering among three parti- 
cles is independent of the order in which the particles 
scatter among one another. For an SU(2) invariant sys- 
tem, the Yang-Baxter relationship reduces to the simple 
condition^ 

b(p,q)/c(p,q)=g(p)-g(q), (2.13) 
where g{p) is some function. In case i), g(p) is given by 
(p-e da - U aa /2) 2 ^ 



g(p) = 2r v 

^± ct^ eta. 

while in case ii), g(p) equals 

(p-e d -U/2f 



g(p) 



(2.14) 



(2.15) 



In the second case, g{p) is equivalent to a single dot with 
level broadening T = J2 a ^ a - 

Having determine under what conditions the dot-lead 
Hamiltonian is exactly solvable, we are now in a position 
to construct A^-particle eigenstates in a controlled fash- 
ion. An eigenfunction with spin 2S Z = N — 2M, is char- 
acterized by a sea of N electrons each carrying momenta 
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{Qi}iLi an d so total energy E = J2i Qi- I 11 a periodic sys- 
tem of length L, integrability allows us to write down in a 
compact form the fc^-quantization conditions (the Bcthe 
ansatz equations): 

AI 

iqjL+iSiqj) = TT TO) ~ X a + V 2 , 

n r X a -g(q j )+i/2 yr X a - Xp + i 

=i X a -g(q,)- l /2 - n Xa -X p - i - [ ■ ) 

These equations are identical to those for the ordinary 
Anderson models but for the form of 5(q) (here given in 
Eqn. |2.4p ). The M-A a appearing in the above equations 
are indicative of the spin degrees of freedom. 

The Bethe ansatz equations are derived in the same 



fashion (using the algebraic Bethe ansatz) as described 
in Ref. fHf). However it is important to also explicitly 
construct the N-particlc wavefunctions. There is an a 
priori possibility that even though the integrability con- 
ditions (i.e. Eqn. (|2.13p ) are satisfied, the N-particle 
wavefunctions posited by the Bethe ansatz are not actual 
cigenfunctions of the Hamiltonian. This possibility has 
been seen previously in Hubbard models with more than 
two species of electrons^ In such a case what spoils in- 
tegrability is the possibility that more than two electrons 
sit on the same site. A similar possibility thus exists here 
where more than two electrons can sit on the dots. 

For simplicity, we only consider the case where two 
dots are involved. The N-particle wavefunctions with M 
i electrons (assuming N > M > 2) in this case take the 
form 



\ipN,M) = I dxi- ■ ■ dx N g(xi, . . . ,xn)c\(xi) ■ ■ ■ c\(x m )c\(xm+i) ■ ■ ■ c\(x N ) 

/~] / dx\ ■ ■ ■ dxM- Sl {o)dxM+i ■ ■ ■ dxN- Sl ( (T )e aa (xi, . . . ,in- Sl ( r ),ijn+i, ■ ■ ■ ,a;jv-s T (cr)) 

a=l,2,<T=Tii 

xcj(a:i) ■ • ■c\{x M _ Sl ( a) )c\{x M +i) ■ ■ ■cj(x A r_ Sr((T) )(i t (TQ 

+ ^2 f ■ " •d®M-s l (a)-8 l (<T')d x M+l ■ ■ •dXlf-s i ,(,<T)-s l (a') 
crcr' aa' 

x /o-o-'aa' (^lj • ■ ■ j XM-s l {a)-s l (cr')7 XM+1, ■ ■ ■ , X N-s r (cr)- Sf (cr' )) 

Xcj(xi) • • ■c\{x M _ Sl(a) _ Sl[a ,))c\(x M +i) ■ ■ ■ c\{x N _^ {a) _ Sl{al) )d) aa d) a , a < 

+ J dXi---dx M -s l {o)-s l {cj')-s l (cT'')dXM+l---dx N _ sy ( rT )_ s ^( rT ^_ s ^( a ii} 

a a' a" eta' a." 

r crcr'cr"aa'a" [Xy, . . . , Xfr[_ s ^ ( ff )_ Si (cr') — S| (cr") ) X M+1 1 ■ ■ ■ > %N— s T (cr)— s T (cr') — Sf(cr")) 

Xcj(.Ti) • • ■c\(XM-s l (a)-s L ((T')-s l (a"))c\(xM+l) ' ' ' c{ (x N - Sl ( CT )-s T (<r')-s T (cr") )d f a a d\> a > d) a " a" 

+ J dxi ■ ■ ■ dx M ~2dx M+1 ■ ■ -dx N -2 h(xi, . . . ,XM-2,XM+1, ■ ■ ■ ,%N-2) 

Xcj(xi) • • ■ c\(xm-2)c\{x m +i) ■ ■ ■ c\-(x N ^ 2 )d^id\ 1 d^ 2 d ii 12, (2.17) 
I 

where = S|(T) = 1 an d s l(T) = s t(1) = 0- The are single particle energies. Under the Bethe ansatz the 

different terms of the wavefunction correspond to differ- coefficients, g{xi), e aa {xi), fa-a-'aa'{xi), v a cr'cr"aoi'a"{xi), 

ent numbers of electrons sitting on the dots. We look for and h(xi) take the form 
eigenstates with total energy E = X^iLi 1i where the q t 



g(x 1 ,...,x N ) = ^ TT 9g„ m (xj) ■ A NtM (xi\o-); 

M-Si(cr) N-s y (cr) 

e aa (x 1 ,...,x M -s l (a),XM+i,---,x N _ Sl r (T )) = ^ 9q aW (xi) g q<j{z) (x l ) 

cr£S N i=l i = M+l 
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M 

n 



N 

n 



°q„(i)B aa N,M{xl, • • • l ^M-SifffJ^AZ+l: • • • j x N-sj(a)W)', 



i=M — Si (ct)+1 i=Af-s T (cr) + l 

facr'aa' {x\ , ■ ■ ■ , X} i ,[-s l (a)-s l (a')i X M+1: ■ ■ ■ , I£jV— s T (cr) — s T (er') ) = 
Af— (it) — si (a 1 ) N— s T (tr) — s T (er') 

x n flfcw^*) n Sflcrw^i) 

o-eSjv »=1 »=Af+l 

M iV 

x n e< ?-w n e< ?-(o 

i=Af— sj.(o-)— S|(ct')+1 i=JV— s-j (<r)— S| (<r')+l 

XCto'bc'^M^Ii ■ ■ ■ , XM- Sl (a)- Sl (cy'),XM+l, ■ ■ ■ , ^Af-s T (<x)-s T (cr') l^); 

r mrV"iia'(i" (cCi , . . . , X^/_ s ^ ^ CT )_ S ^ ) — (<r" ) ! ^-M+lj • • • i X N—Sf (tr) — s T (<r') — s 1 (""") ) = 

A/-s i (cr)-s i (cr')-s i (cr") iV-Sf (<t)-S T (<t') -Sf (<r") 



i=M+l 



z— Af— S| (cr) — sj.(c') — sj (a -/ ')+l i—N—Sf (a) — s-[((T') — sy(a") + l 

XD (T(r i (r n aa t a "N,M(xi, ■ ■ ■ , aJM-8j.(o-)-s J .((7 / )-8j.(c7")> a; Af+li • ■ • ) a; iV-s T (cr)-s T (cr')-ST(' :r ")l Cr )' 

M-2 JV-2 
,Xm-2,2;M+1, • • • ,ZN-2) = ^2 e< ?-W II ff<M<) ( x i) 



aeS N i=l i=M+l 

xEn,m(%i, ■ ■ ■ ,X M ~2,XM+1, ■ ■ ■ ,XN-2\cr). 

I 



(2.18) 



Here ^o-eSw * s a sum over the permutations of N. g q (x) 
is the electron portion of the single particle eigenstate de- 
fined in (|2.3[) with energy q. Similarly e q is the dot por- 
tion of this same wavefunction. The coefficients -Ajv.ju 
through Ej^ m depend upon both a and the relative or- 
dering of the Xi (for example, An,m only changes value 



B. Analysis of Bethe Ansatz Equations 

1. Description at T — 



N-particle states with spin projection, 2S Z = N — 2M, 



as the Xi cross one another). Their forms are analogous of the dot-lead system in general consists of a multitude 



to those of the Hubbard model detailed in Rcf. 0) 



By direct substitution we check that this posited form 
of the wavefunction does indeed correspond to an eigen- 
function of the Hamiltonian. While straightforward, the 
associated algebra is tedious and so is suppressed. 



of different types of solutions for the q's appearing in the 
Bethe ansatz quantization conditions of Eqn. (|2.15p . But 
at T = only two solutions are relevant for the ground 
state. The ground state at T = is composed of 1) N-2M 
real values of q and 2) 2M complex values of q interlinked 
to M real values of A according to the rule 

q ± (X)=x(X)±iy(X); g(q± (A)) = A ± i/2. (2.19) 

In the continuum limit, densities p(q) and er(A) describing 
the distributions of q and A respectively may be derived^ 



While we have only consider two dots in parallel, we 
can easily generalize to N-dots. In order for N-dots to be 
integrable we again require that V\ a /V2a = Via' /^2a' for 
all a and a 1 and that constraints analogous to those in 
Eqns. (|2.7p and (|2.8[) are satisfied. This type of gener- 
alization is akin to the multi-lead generalizations of the 
single level Anderson model found in Refs. [45| and[46l 



P(q) 
a(X) 



^ + ^ + 9'{q) [ dX ai (g(q) - A)<r(A); 
2tt L J q 



a/ (A) A(A) 



TC 
B 



L 



Q 



dX' a 2 {X' - X)a{X') 



dqai{X - g{q))p{q), 



(2.20) 



D 
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where 



Mm) 

A(A) 
a n (x) 



^-d q 6(q); 
Zir 



d x ReS(x(X)+iy(X)); 

ir 

^-d x 9 n {x) = — 1 . 
2ir ir (n 2 + 4x 2 ) 



(2.21) 



In the equations for the charge, p(q), and spin, cr(A), 
distributions appear the Fermi surfaces, Q and B, and 
"band" bottoms Q and —D. These mark out the range 
of q and A over which excitations appear in the ground 
state. — D is the lower band edge of the charge exci- 
tations. As each A has two associated complex q's, we 
expect q q+ (Q) + q q _(Q) = 2x(Q) = —2D, thus deter- 
mining Q.— The Fermi surfaces Q and B are fixed by 
insisting that the overall spin and particle number, 



N - 2M 



L I dqp(q); 

l-D 



M = L 



d\a(X), 



(2.22) 



Q 



are reproduced. 

In understanding transport properties we will want to 
construct electronic-like excitations (even if zero-energy) 
about this ground state. These excitations are in turn 
related to specific solutions, q and A, of the Bethe ansatz 
equations. If for example we add an electron to the sys- 
tem, i.e. take N — > N+l, to effect this change, we must 
add a real g-state as well as a A-holc. Thus we must 
parameterize the energies, e q (q) and e(A) of the these ex- 
citations. Equation s g overning these energies are easily 
derivable (see Ref. (|36T and are given by 



hierarchy of excitations of the Bethe ansatz equations 
must be considered. The two coupled integral equations 
at T = are replaced by an a correspondingly infinite set. 
These equations are closely related to those analyzed in 
Ref. (HfHII). However some important differences are 
present and so we give the description of the system at 
finite temperature in some detail. 

At finite temperature, the possible excitations of the 
system fall into three classes 

i) real q: As we have already seen, such solutions to the 
Bethe ansatz equations appear in the zero temperature 
ground state at finite H. 

ii) n-spin complex with 2n associated complex q's: 

The n-spin complex groups together n different A's ac- 
cording to the rule 



X nj = y + 



1 



- j), 3 = !>••• ,n. 



Here A™, the centre of the complex, is a real number. For 
each A™ 5 , there are two associated q's: 

9{q +n3 ) = A" + i(| + l-j); 
g(q-^) = A"+i(|-j). 

The simplest of this type of complex (n — 1), as we have 
already seen, fills the H = T = ground state. 

iii) n-spin complex with no associated q's: As in 

case ii), the spin complex groups together n different A's 
according to 



j), J I- 



. n. 



ca(A) 



H 



dXe x (X)a 1 (X - g(q)); 



2x(A)- / dA'e A (A / )a 2 (A / - A) 

B 

+ 1 g'(q)e q (q)a 1 (g(q)-X). (2.23) 

J-D 



e q (q) is a monotonically increasing function that is equal 
to zero at the q- Fermi surface, i.e. e q {B) = 0. For q > B, 
e q (q) > indicating it costs energy, e q (q), to add a state 
above the q— Fermi surface while for q < B, e q (q) < 0, 
indicating it costs energy to add a hole below the Fermi 
surface. Similar consideration hold for e\ (although here 
e(A) is a monotonically decreasing function which passes 
through at A — Q). 



2. Description at Finite T 

The description of the system at finite temperature 
is considerably more involved. At finite T, an infinite 



However here there are no associated q's. 

These excitations are governed by the following set of 
density equations describing their occupancies at finite T 
(akin to Eqn. (|2.16p for the simpler T = case) 

p P (q) + Ph(q) = ^ + 

+ 9'(q) / dXa n(g{q) - A)(<r pn (A) + <y' pn {X))- 

n=1 J-oo 

m <(A) A n (A) 

07m (A) = 1 

TT L 

a n (X- g(q))p p (q) 
~ / dX'A nm (X - A')crpm(A'); 

m=1 J-oo 

oin(A) = / an{X- g(q))p p (q) 



— OO 

oo 



E 



dX'A nm (X — X')a'(X'); (2.24) 
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where x n (X) and A„(A) are denned by 



, U 



x n (X) = VWTRe (X + i-) 1 / 2 +n{- + e d ) 
An (A) = --d x 5 n (X) 



—9a Re S(—VWT(X + i^) 1 ' 2 + U/2 + e d ) 



n-l 



k=l 



T-dx V 5(-VzUT(X + Un - 2k)) 1 / 2 + U/2 + e d ) 

Z7T *■ — ' 1 



+5{V2UT (A + - [n - 2k)) 1 / 2 + U/2 + e d ))j, (2.25) 

The kernels in the above density integral equations are 
2n 1 



Q"n(X) — 

it [n l + 4A Z ) 
A nm (X) = 5 nm S(X) + a|„- m |(A) 

min(n,m) — 1 

+2 a \ 

7i— m\ -\-2k (A) + (A). (2.26) 

k=l 

In the above, p p /h(q) describe the particle/hole occu- 
pancy of the real q— excitations. The total density of 
states, p(q), is equal to the sum of the two 

p(q) = Pp(q) + Ph{q)- 

To compare with the T = density equations, we note 
that at T = 0, 

p p (q) = e(B-q)p(q); 
Ph(q) = Q(q-B)p(q). 

&p/hn(X) similarly describe the finite temperature parti- 
cle/hole occupancies of the n-spin complexes with associ- 
ated complex g's (here A refers to the string centre). At 
T = we have for comparison 

<7pn{X) = 0,71 > 1; 



0pi(A) = f7i(A)6(A-Q); 
ofci(A) = ( 7 1 (A)9(Q-A). 



(2.27) 



Finally v' p p ln {X) describes the particle/hole densities of 
the n-spin complexes without associated complex-q's. At 
T = 0, <j' pn {X) are uniformly zero (none of these excita- 
tions appear in the ground state). 

The energies of creating the various types of excita- 
tions are indifferent to the form of the impurity scattering 
phase (as at T = in Eqn. (|2.2ip ). Thus these energies 
are the same as for the single dot case (see Ref. (|36T)) 
and are given by 

e q(q) = 1 



tJ2 [°° dXiog S e ^ K(A-g(g)); 

^[J-oo /(-£An(A)) 



log(/(e An (A))) = --x n (X) 



m— 1 

logCMJA))) 



dqg'{q) \og{f(-e q (q)))a n (g(q) - X) 
J dX'A nm (X - A') log(/(-e Am (A'))); 



dqg'(q) log{f{-e q (q)))a n (g(q) - A) 

+ J2 f dX'A nm (X - A') log(/(-e Am (A'))), (2.28) 

where /(e) = (1 + cxp(e/T)) _1 is the Fermi distribution. 
These energies are related to the various particle-hole 
densities via 

P P /h{q) = (Pp(q)+ Ph(q))f(±e q {q)) = p(q)f(±e(q))\ 
&p/hn (A) — (o>i(A) + a hn (X))f(±e\ n (X)) 

= <r(A)/(±e An (A)); 
a p/hn( x ) = Kn(A) + cr' hn (X))f(±e Xn (X)) 

= a'(X)f(±e' Xn (X)); (2.29) 

These equations are arrived at from minimizing the free 
energy of the system. The solutions to Eqns. (|2.24[) and 
(|2.28p represent the distribution of excitations in thermal 
equilibrium. 



III. COMPUTING TRANSPORT PROPERTIES: 
METHODOLOGY 

Here we describe how to extract the scattering phases 
of electrons off the dot system from the Bethe ansatz 
equations. We do so both for the system at zero and finite 
temperature in equilibrium. We believe that the results 
are obtained by this methodology are exact at T = and 
an excellent approximation at finite T . The nature of this 
finite T approximation is discussed in detail in Ref. (f36h . 



A. General Considerations 

To compute transport properties we must translate the 
effects of the map between the L/R and the even/odd 
electrons. Excitations in the even/odd picture scatter off 
the collection of dots with some pure scattering phase, 
5 e / (E), where E is the energy of the excitation. This 
is true regardless of energy as there are no inelastic pro- 
cesses to take into account. The excitations are exact 
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cigcnstatcs of the system and not merely asymptotic ap- 
proximations. The transmission/reflection (T/R) ampli- 
tudes of an excitation in the L/ R basis are related to the 
two scattering phases, S e / via the following ansatz: 



iS e (E) 



K(E)+T(E); 
1 = K(E) -T{E). 



(3.1) 



5 (E) is necessarily as the odd sector is uncoupled from 
the dots. 

The computation of the 5 e (E) for a given excitation 
is described in detail in Ref. (|36l). It is based on a 
technique used to compute the magnetoresistance in the 
Kondo model42 The scattering phase is identified with 
the portion of an excitation's momentum scaling as l/L: 



p = p hulk + p lmp /L = p bulk + 5 e {E)/L 



bulk 



(3.2) 



For the purposes of computing transport properties, we 
are interested in the scattering of electronic excitations. 
In the Bethe ansatz solution of this model, all excita- 
tions see a spin-charge separation. In order to construct 
an electronic excitation, we must glue together a charge 
excitation with a spin excitation. How this is done and 
potential (finite energy /temperature) pitfalls are again 
discussed in detail in Ref. (|36l ). A spin f electron can be 
treated as a composite of a right-moving real g-excitation 
and a left moving A-hole (here by a A-hole we mean a hole 
in a 1-spin complex with associated complex q's). Thus 

4lectronT=Pg mP (9)+f I A nlP W- (3-3) 

The energy of this composite particle is 

e c i = e q (q)-e x (X). (3.4) 

The scattering of spin J, electrons can be accessed through 
a particle-hole transformation. 



B. Zero Temperature 

At T = 0, we are interested in zero energy excitations 
at the Fermi surface. The conductance through the dot 
is given by: 



G = -(\T t (E = 0)\ 2 + \T l (E 



5 el (E = 0) 



2 



) + sin 2 ( 



))• (3-5) 



The zero energy composite excitation is constructed by 
choosing the q— , A-cxcitations at their Fermi surfaces q = 
B/X = Q. (At zero temperature e q (q = B) = e\{X = 

Q) = o.) 

A simplification in the computation of scattering 
phases arises from a relationship between the impurity 



momenta and the impurity densities. The impurity den- 
sity is simply the portion of the density scaling as l/L, 
that is, the piece sensitive to the presence of the dots: 



p(q) = Pb u^(q)+P imp (q); 
<t(\) = a bulk (A)+a imp (A). 
Thus p imp {q) and cr im P(A) satisfy 



(3.6) 



p imp (q) = A(q)+g'(q) / dX ai (g(q) - A)a im P(A); 



(3.7) 



Q 

a imp (A) = A(A) - / dX'a 2 {X' - A)a lmp (A') 

,-B 



-D 



dq ai {X-g(q))r P (q), 



In Ref. (|36f ). it is then shown that 

d qPq mp {q) = 2V mp (?); d xP T P W = -2™ im P(A). 

(3.8) 

Thus by computing the impurity densities we can access 
the scattering phases of the excitations. In particular we 
see 



(3.9) 



p™P(B) = 27T / dqp^(q) 



-D 



p™ p (Q) = 2tt / dXa lmp (X). 



We point out that the integrated densities do not neces- 
sarily equal the number of electrons sitting on the dots. 
For example 

T^i + Hi^ l B dqp imp (q)+ [ Q dXa imp (X), 

I J-D JQ 

as would be the case for a single dot. Rather the r.h.s. of 
the above equation includes contributions coming from 
l/L changes in the electron density of the leads. We will 
discuss this further in Sections V and VII. 



C. Finite Temperature 

At finite temperature, the calculation is less straight- 
forward. We must both take into account all the ther- 
mally excited excitations and their dressing of the elec- 
tron scattering phased For a variety of technical reasons, 
this is only feasible at the symmetric point of the dots, 
i.e. Uu/2 + €di = as explained in Ref. ((361 ). At finite T 
the electronic excitations have a Fermi distribution and 
so the conductance, G, equals 

G = ^ J d E(- ^)(\T t (E)\ 2 + \T l (E)\ 2 ). (3.10) 
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To determine the scattering probabilities, Tj n (E) = 

sin 2 ( ^j-^- ) , we construct excitations above the thermal 
ground state discussed in Section II. B. 2 in analogy to 
what we did at T = 0. The electronic excitations are 
built out of real q and an n = 1 spin-charge complex, A 
(in the language of Section II. B. 2). We have specifically 



(3.11) 



However here p]^ np (q) and p I 1 mp (A) are dressed by the 
presence of the other thermal excitations in the ground 
state. The finite temperature form of p lrap (q) andp I 1 mp (A) 



,.36 



°° />oo 

p^(q) = 5(q) + Y / / e n (g(q)-X) 

1 J — OO 

x«P(A)+4n mp (A)); 



71=1 • 



p£ p (A) = 26 1 (\)+ / dqp^(q)0 1 (X-g(q)) 



OO „ 

£ J dA'E lm (A-A')<P(A'), (3.12) 



where the kernels n (X) and <7i m (A) are given by 

6JX) = 2tan _1 (-A) -it; 

n 

S, mi (A) = S| n _ TO |(A) 

min(n,m) — 1 

+ 2 X] ^|n-m|+2fe(A) + 0„+ m (A). (3.13) 
fc=l 



Comparing the above with Eqn. (|2.25|) , we find the rela- 
tionship, 

d q p i ^(q,T) = 2np imp (q,T); 

d x p^(X,T) = -2na i T P (q,T) (3.14) 
Hence at finite temperature 

Sel^T) = 2tt r dqp™*{q,T) 

J-D 

rQ 

+2?r J dAal mp (A,T), (3.15) 



a generalization of Eqns. (|3.3p and (|3.9p . 

We need to determine how to choose q and A at finite 
temperature. If we are interested in constructing an ex- 
citation of energy, e e i, out of a q— particle and a A-hole 
we must insist that 



f-el = £q{q) - e(A). 



(3.16) 



This would seem to provide a two parameter space (q, A) 
for characterizing excitations. However as explained in 
Ref. HH), we fix A to some value Aq allowing only q to 



vary. One justification for doing so lies in the nature of 
the physics at the symmetric point. At the symmetric 
point, Kondo physics will be present in some form gov- 
erned by a scale, Tk- Of the two energies only e(q) sees 
variations on this scale. e(A) in contrast see variations on 
the much larger energy scale, y/UuTi governing charge 
fluctuations. Thus it is natural to only vary q. To fix 
Ao we exploit the fact that at the symmetric point of the 
dot system, we expect the scattering phase to similarly 
symmetric about zero energy, i.e. 



5 e i(e,T) = S e i(-e,T). 



(3.17) 



We thus have 



S el] (e eh T) = 2ir / dqp lmp (q,T) 

J-D 



-2tt / dAal mp (A,T), (3.18) 

'An 



with q chosen such than 

£ei = e(q) - e(A ). 



(3.19) 



To compute <5 e f from Eqn. (|3.18[) . we need to compute 
p lmp (q) and er" np (A). This can by done by solving Eqns. 
(|2.24p and (|2.28p . However at the symmetric point and 
at energies much smaller than y/UuTi, these equations 
can be simplified. If we recast the energies as 



0„(A) 



-e n (X- -log(— )); 

1 TT 1 

-^e(-5(e) + -iogt)); 

1 7T 1 

^4(-A + -log(— )); 

1 7T 1 



A = 



V2UT 



(3.20) 



we can show that they satisfy the following set of integral 
equations 



&(A) = -/ dX's(X- A')log(/(TC„-i(A))/(Te i+ i(A))) 

J — OO 

-5 nl e* x , (3.21) 

where = <j> n or <f>' n and s(A) = cosh - 1 (tt A) /2. 

The impurity density equations take the initial form 
given in Eqn. (|2.24|) . If we focus upon energies <C 
y/UuTi, we can safely assume 

47(A) = 0; ^7 = 0; m > 1. (3.22) 

The impurity densities then reduce to 

p imp (q) = A(q) 

+9'(l) £ f°° dXa ^) - A)< np (A); 

m 1 J — OO 
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+g'(q) / d\ai(g(q) - X)a^ p (X) 



<7(A) = Ai(A) - / dqp™ p (q) ai (X - g(q)) 

J — OO 

d\'a 2 (\-\')a™ p (\'); 
^T P (A) = / d qP f^( q )a n (\-g(q)) 

J — oc 

/oo °° 
dA' J] A„ m (A - A')^p(A'). (3.23) 
-°° rn=l 



Using the inverse of the matrix A nm , 

4m(^) = 6 nm S(X) - s(X)(S nm+1 + <W_i), (3.24) 
together with the relationship 

/oo 
dA'A^(A- A> m (A' - A"), 
-OO 

(3.25) 

we obtain 

/oo 
dAs(A- 5 (g))A 1 (A) 
-OO 

/>co 

/ dA,s(A-. 9 ( g )K7 P (A); 



%7(A) + *Z P W = d\'R(X A')Ai(A') 



-Ai(A') 



d9/?r p (9)s(A -#(?)); 



<r(A) + ^T P (A) = «nl / dqp™ p (q)s(X - g(q)) 



+ d\'8(\-\')(oZZiW + v%£M) (3-26) 



In this form, the equations are easily solved numerically 
through iteration. 

This computation of the finite temperature conduc- 
tance comes with some potential pitfalls as discussed in 
Ref. (f36h . However as shown in this same reference we 
are able to accurately reproduce the universal finite tem- 
perature linear response conductance scaling curve for a 
single dot. Uncertainties arose with the application of 
this methodology because of uncertainties with identify- 
ing the correct bulk scattering states. Here however the 
choice of bulk scattering states is exactly the same - dif- 
ferences between the treatment here and Ref. (f3tl ) only 
arise at the level of impurity scattering. As our choice of 
bulk scattering states was unproblcmatic previously, we 
expect it to be similarly so here. 



IV. ANALYSIS OF EQUATIONS GOVERNING 
ZERO TEMPERATURE TRANSPORT OF 
DOUBLE DOTS IN PARALLEL 

In this section we will analyze in detail the equations 
governing zero temperature transport in double quantum 
dots. As discussed in Section II. B, there are two sets of 
constraints under which a pair of dots in parallel is inte- 
grable. The most interesting set of constraints is given 
in Eqn. 



U a + 2eda — U a > + 2ed a >; 



(4.1) 



These conditions correspond to well separated parallel 
dots. They are neither capacitively nor tunneled cou- 
pled. The second set of conditions, corresponding to the 
constraints, 



II- T — F /• 



(4.2) 



leads to transport that is identical in nature to a single 
level dot. As such we will focus primarily in this section 
on parallel dots meeting the first set of constraints and 
only briefly discuss dots constrained by the second set at 
the section's end. 

At zero temperature, as discussed in Section III. A, the 
scattering phases which determines the transport are ex- 
pressed in terms of the two impurity densities, p lmp (<7) 
and er lmp (A). These two densities satisfy the two integral 
equations (Eqn. (3.7)): 



P imp (9) 
cr imp (A) 



Po mP (9) +9'(q) d\ ai (g(q) - A)a lmp (A); 



Q 



Q 



a mp (X)- dX'a 2 {X' -X)a lmp (X') 
JQ 

dqa^X-g^p^iq). (4.3) 

'-D 

The 'Fermi surfaces' of the dots, Q and B are determined 
by the conditions given in Eqn. (2.19). As Q and B (up 
to 1/L corrections) are given solely by the bulk densi- 
ties <7 bulk (A) and p hulk {q), their behaviour, as a function 
of the two quantities 2edi/2 — an d r^/Tn^, is 

exactly the same as that of a single level dot and so dis- 
cussed extensively in Refs.dH) and ((361). 

The equations in Eqn. (|4.3p can only be solved in gen- 
eral numerically. However they do admit analytic solu- 
tions in certain cases in two limits: a) |e<ji — e^l 3> Ti, T 2 
and b) \edi — e^l <C ri,r 2 . In both cases the source 
terms, p 0lrap (q) and <7 0lmp (A) can be written in the form 



imp / \ 
imp / \ \ 

c (A) 



AoT P (9) 



ATO; 



where 



A ? p (<?) 



'01 



r, 



'(A) 



'02 



'(A), 



(4.4) 



f? + (^e dl ) 2 ' 
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a™ p (X) = -Be-d x A 0i (x(\)+iy(\)y, 



dqA 0i (q)ai(X - g{q)) 

1 

-ai(A - g(e di + if;)) + (3 l a 1 ([3 l X - 7.;); 



Ui 
2 

2a,T,; 



oLi — — + tdi — tdi ; 



Pi 



(4.5) 



The parameters Idi and I\ are case dependent. For case 
a), we obtain: 



kdi - £d 2 | > r 1; r 2 : 

Tl/2 = I\/2 ' 

£dl/2 = e dl/2 

and for case b) 

kdi - £d2| < ri,r 2 : 



rir 2 



(edi - £d2) 2 

rir 2 

e dl/2 — £d2/l 



(T1/2 - r 2 /i); 



(4.6) 



ri = r! + r 2 ; 

Cdl + £ci2 
e dl — 7: • 



p (edi — ed2) 2 rir 2 
2 ~ (r! + r 2 ) 3 : 

£d2 = £d2 + (fdl — £d2) 



ri + r 2 



(4.7) 



If £di = £d2 exactly, then r 2 vanishes and Ao 2 (g) is iden- 
tically zero. 

Given the division of the source terms into two pieces 
and the linearity of the equations in (|4.3|) . we can cor- 
respondingly write er(A) = ci(A) + (72(A) and p(X) = 
Pi(X) + p 2 (A) where 

fT'(q) = pZ p (q)+9'(D [ Q d\ ai (g( q )-\yrW; 

ar P (X) = a^(X)- [ Q dX'a2(X'-Xyr(X / ) 



dqai(X- g{q))p i p (q) 



(4.8) 



A. Linear Response Conductance at Zero Magnetic 

Field 



We compute the Fermi surface scattering phase, 
5 e , governing the linear response conductance, G = 
2e 2 //isin 2 (<5 e /2), via 



S R = 2tt / K mp + ( r 2 mp )(A) 

Jq 



J — 00 

2t7 - S el - 5 e2 . 



imp 



)(A) 



(4.9) 



The quantity, , in first line of Eqn. (|4.9p is the impu- 
rity particle number while the corresponding integrated 
quantity in the second line is the impurity hole number. 
cr™ p can be obtained via a Weiner-Hopf analysis of Eqn. 
(4.8) (with B = —D). Our separation of the impurity 
densities into two pieces allows the computation to pro- 
ceed along the lines of Ref . (f35h . We relegate the anal- 
ysis to Appendix A simply reporting the results for 8 e i 
below. These results apply to both 8 e \ and S e 2 provided 
f-di =/= (-d2- If instead td\ = £d2, $e2 is identically zero and 
the equations below apply only to 5 e \. 

The analysis of 5 e i breaks down into three cases de- 
pending on the value of Q, the parameter describing the 
position of the Fermi surface of the A-excitations. 

Case i: Dots Near the Particle-Hole Symmetric 
Point Q < 0: Near the symmetric point, i.e. U{/2 + 
e d i < VUT, then Q < 0. Setting J, = I7 l - Q, where 



17 1 = 



2U,T,, 



5 e i is then given by 



5 ei = Kn + K i: 



(4.10) 



(4.11) 



Here Kij is the contribution to the scattering phase, d e i, 
coming from the bare impurity densities, A™ p ((/) and 
(Tq™ p (A). In this case we find (see Appendix A) 

(-1)™ (n+l/2\ 71+1/2 1 



' n + 1 



n=0 



/2 



r(n + 1) 



x / dqA m { q )e-* {2n+1){a{q) - Q) ; 

J —OO 

K i2 = Bign(l-A)V5rP 



°°dujfujV e - 2 ™ J < 



o lu \e J r(l/2 + w) 
sin(27rw(7 4 - 1/2)) 
003(7101) 

+ sign(l-A)V^> v 

S) n + 1 / 2 V e / r(n+l) 

xcos(7r(7 < -l/2)(2n + l)), (4.12) 

where P indicates that the principal value of the integral 
is to be taken and 7, = 1/2(1 — 07 )sign(l — /?;). The 
value of the Fermi surface, Q, was found in Refs. |35h 
and (|36h to be determined implicitly by 



TTvUs G (z7r(2n+ 1)). 



X/2TAT7 



7T — (2n+l) 3 / 2 

n=0 V ; 



(4.13) 
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Case ii: Mixed Valence Regime of the Dots: Q > 

and J,* = /r 1 — Q > 0: For these values of the parame- 
ters, we are in the mixed valence regime of the dots. Here 
8 el is given by 



K a + K i2 



K a = 7r(V2-l-^J l + n^J? 



+ A^iPi 1 -1)4 + 0(J?) + 0(MP^ - i))^ : 



K l2 = sign(l-A)V^P - - 



/ uj \ e J r(l/2 + w) 
sin(27rcj(7i - 1/2)) 
cos(7ro;) 

+ sign 1 - Pi)V*2l ^TTT^ r7 ZT> 

^ n + 1 / 2 \ e y r(n + l) 

xcos( 7 r(7 i -l/2)(2n+l)). (4.14) 
If Ji is small, if,2 can be approximated by 

K i2 = -sign(l-A)— tan" 1 (4.15) 

For 74 small (as would be the case if |e<ji — ed2| IY2) we 
see that -Ki2 undergoes a rapid variation as Jj is varied. 
This in turn will imply a rapid variation in both the T = 
linear response conductance and the number of localized 
electrons in the system as the gate voltage is swept. 

In this case Q needs to be determined by a self- 
consistent numerical evaluation of Eqns. (|2.23p . If Q ~ 1 
however, it is a reasonable approximation to take 



Q = q* + — ln(27reg*), 
where q* = (e di + C/,/2) 2 /{2U l T i ). 



(4.16) 



Case iii: Empty orbital regime: Q > and Jj = 

J. — Q < 0. In this parameter range, the dot filling is 
small. S e i is then given by 

Sei = Kn + K i2 

Ka = ,4r^r(l/2 + W )W>- 

x sin(7rw(l + — )); 

Pi 

1 POO 1 

K l2 = 7rsign(l — Pi) H — p / — r(l/2 + o;) 

V 71 " w 



i 2 ™ Ji sin(7ra;7 i ). 



(4.17) 



If Ji is small, the expression for K a in Eqn. (|4.15p holds 
equally well. Similarly, in this regime, Q is given as above 
in Eqn. (j4~TrJ)l . 



B. Finite Field Transport 

In finite field, we must compute the scattering phases, 
fief/ii separately for spin up/down electrons. Like at 
H = 0, these scattering phases are given in terms of 
the impurity densities: 



s e] = 2^ / dq^r^+pni)) 

-D 



+2W (ar+^ mP )(A); 
JQ 



S el = 27T / (cr™ p + (T™ lp )(X). 



(4.18) 



In general these equations cannot be solved analyti- 
cally. However at the symmetric point they do admit 
such a solution. At this point ^ Q = 00 and Oi can be 
expressed solely in terms of pi (simply by inverting Eqn. 
4.8 via Fourier transform): 

/oo 
dA'(l + a 2 )- 1 (A-A'K7(A') 

dqf/r P (q)s(\-g(q)). (4.19) 

1 

In turn, we obtain an equation for /9™ p independent of 



Pr P (l)=Pr(l)-9(q) / dq'R{g(q)-g{q'))p^{q') 1 

J — OO 

(4.20) 

where 



1 f 

and we can write p™ p (q) as 



(4.21) 



p% p (q) = Aoi(q)+g'(q) / dqR(g(q) - g(q'))A oi (q) 

J — OO 

/OO J 

x —-r . (4.22) 

2cosh(f) V ' 

We are most interested in the behaviour of p - mp for q < 
as this is the region relevant for magnetotransport when 
H is on the order of the Kondo temperature. Following 
Ref . HH) , for q < the first two terms may be rewritten 
as 

/OO 
dqR(g(q) - g(q'))A 0i (q) = 
-OO 

r°° i 1 ,e- ,u (9M- f . :1 ) + ¥ 1 -r 1 

-g'(g) =■- j—TTzr, • (4.23) 

/.oo 2tt 2cosh(2-) 
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At the symmetric point, the scattering phases similarly 
simplify reducing to 



J-D 

S el = 27T-n [ B dq(p i r(q)+ P T P (q)). (4.24) 
J-D 

We thus see that even though <5 e j differs from 8 e i, the 
conductance for each spin species is the same. 

With the source term in this form, Eqn. (|4.23[) can 
be solved by a Wiener-Hopf analysis. We again relegate 
the details to Appendix A and here simply report the 
results. Firstly we consider the dependence of the limit 
B, the Fermi surface of the q-excitations, on the magnetic 
field. It has the same functional dependence on Ui + 2tdi 
and TiUi as in the case of a single level dot. Thus we can 
simply borrow results from Refs. (f35h and (HF 



H 

2^ 



8tt 2 



1/2 



duit 



1 



(ioj + S) 1 / 2 uj - iS 



±) V 2?re 



(4.25) 



In the small H limit, one can deform the contour into 
the lower half plane, merely taking into account the pole 
nearest zero. We so obtain 



B 2 1 , fireH 2 
= log 

2U l T l 2tt 6 V 41T 



(4.26) 



Next we write down expressions for J_ D dqp 1 ^ np (q) 



dqpT P {<l) = ^f 0nd ° + ^ CharS °- 



(4.27) 



D 



Here S' l Kondo encodes information that varies on the 
Kondo scale. It is given by 



S 



20(A - 1 



2n+l , x n+i 

2±l) cc . s(7r(n+ i )(!_£-!)), if b >T ~\ 



r(l+n)(n+i) ^t£ kk1 '2 3 / 2 

1 ~ fo°° 77 ( B{T LP 2 ) + w) sm(^) cos(^(l - /37 1 )) 



(4.28) 



if b < ir 1 . 



where the Kondo temperature is defined by 

1 /2 

T* KKY = e~^\ if pi > 1. (4.29) 

(We explain the label 'RKKY' in Section V.) This def- 
inition is well posed as of P± and P2, only one will be 
ever greater than 1. This Kondo temperature is similar 
to what is found in a single level dot but for the form 
of 77 (see the definition in Eqn. (|4. 10|) involving pa- 
rameters renormalized by the effects of interference be- 
tween the two dots) . We can evaluate S^ ondo i n the limit 
H > lf KKY (but H < y/2UiTi/Tre): 



S 



Kondc 



29(ft - 1) 



1 - 



-O 



-O 



2\og(H/T™ KY ) 

/ iogiog(g/r£ KKY u ^(l-ff- 1 ) 2 

I log 2 (H/f£ KKY ) J 8\og 3 (H/f^ KKY ) 
(1-/37 1 ) 2 log log(ff/f « KKY _) \ 



\og A {H/f* KKY ) 
(1-/37 1 ) 4 



\log 5 (H/f^ KY )J 



(4.30) 



where f% KKY = |=t£ kky . We have included the lead- 
ing order term in (1 — P^ 1 ) 2 as this governs the differing 



behaviour of Sf~° ndo in the two cases of je^i — € d2 \ <^ I\,2 
and \e d i - e^l > Fi,2- 

The second term, 5' char g e ; varies on a scale comparable 
to the charging energy of the dot, y/U^l. In the dots' 
Kondo regime this scale will be much larger than Tjc. It 
is given by 



-^charge 



1 \ - e 



^ o r(l + n)(n+i) V e 
x / dqe'^^^+^Ao^iq), (4.31) 



and for fields H ~ yRKKY ma k es a f ar sma ll er contribu- 
tion than that of S Kondo . 



C. Two degenerate dots in parallel 

We now briefly consider two degenerate dots in parallel 
satisfying the conditions given in Eqn. (|4.2[) . This case is 
less interesting than the one described by the conditions 
in Eqn. (|4. 1(1 as it is equivalent to a single level dot 
with Coulomb repulsion, U, and linewidth T = ^2 a T a . 
This follows directly from Eqns. (j2~4"| and (j2~T5")) . the 
equations governing the bare impurity scattering phase, 
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5(q), and the function, g(q), controlling the intcgrability 
of the dot model. Because U aa > = U and = e<j, we can 
make the replacement T — > ^ Q r Q and obtain equations 
corresponding to a single level dot. 

With this equivalence, we can directly apply the anal- 
ysis of the linear response conductance contained in Ref. 
(f36h . And so none of interesting phenomena identified 
in the coming two sections remains. In particular the 
RKKY Kondo effect, the presence of interference phe- 
nomena, and a non-trivial form of the Friedel sum rule 
are all absent. 



V. EXAMPLES OF ZERO TEMPERATURE 
TRANSPORT OF DOUBLE DOTS IN PARALLEL 

In this section we consider various aspects of zero tem- 
perature transport in double quantum dots both in the 
absence and presence of a magnetic field. We will see that 
there is a mixture of different types of Kondo physics with 
interference effects. We will also will sec that the number 
of electrons displaced in attaching the quantum dots to 
the leads can be negative, a contrast to what is found 
in single level quantum dots. We begin by examining 
H = conductance of the dot system as a function of 
gate voltage. 



A. Zero Field Linear Response Conductance as a 
Function of Gate Voltage 

Here we study the behaviour of the T = linear 
response conductance, G, as a function of gate volt- 
age keeping the distance between the dot levels, i.e. 

— £ci2 1 , fixed. Thus we allow the gate voltage to move 
the two dot chemical potentials in unison. In Figure 3 
we plot four examples of G for dot systems with differ- 
ing separations ranging from \e d \ — £<i2\ » Ti,r2 to 
|e<ii — £<12\ = 0. Provided that je^i — e^l ^ 0, the conduc- 
tance exhibits three basic features: i) an RKKY Kondo 
effect at the particle-hole symmetric point of the dot; ii) 
a standard Kondo effect corresponding to a value of gate 
voltage when one electron sits on the dot; and iii) in- 
terference effects whose strength is proportional to the 
separation of the two bare dot levels. 

To understand this better let us begin by consider any 
of first three panels of Figure 3 (a, b, or c). Here |e<a — £d2| 
is nonzero. The linear response conductance trace begins 
at the particle hole symmetric point of the dot where 
(uniformly) e<ji = — 10Fi and e<j2 = — 5Fi, — 3Fi, — 2Fi 
respectively for panels a), b), and c). At this point ex- 
actly two electrons sit on the two dot system. By the 
Friedel sum rule the corresponding scattering phase is 
S e = 2ir and the conductance vanishes. Although the 
conductance vanishes as this point, and counter to the in- 
tuition we have in a single level dot where Kondo physics 
is associated with a unitary maximum in the conduc- 
tance, Kondo physics is still present. We determine its 



presence by studying the low energy density of states of 
the impurity double dot system. We find what amounts 
to an Abrikosov-Suhl resonance. 

In the Bethe ansatz language this amounts to study- 
ing the quantity p lmp (q(e)) where p lmp (q) is the impurity 
density of states of the charge excitations and e q (q) is 
their corresponding bulk excitation energ y 36 i 47 . In zero 
magnetic field, p lmp {q) is equivalent to p 1 ^ np {q) computed 
in Appendix (A3): 



p imp (q < 0) = -2g'(q) 



cos(pn) cosh(n(g(q) — l i 1 )) 
cosh(27r (g(q) - I' 1 )) + cos(2tt/i) 



(5.1) 

where g(q) = {q- e d i - Ui/2) 2 /(2U 1 T l ) and p is given by 



/' 



Kl-^- 1 ) if/3i>l; 

ia-zv 1 ) if/3 2 >i, 



(5.2) 



where (3\^ are defined in Eqn. (|4.5[) . Given the equiva- 
lence of the bulk equations between the single level dot 
and our two dots in parallel, the bulk energy of ch arg e 
excitations, e q (q), can be taken from Refs. (|35[ ) and (|36T ) : 



£9(9) 



_, 

e 



9(g) 



(5.3) 



This relationship is valid provided e q (q) is on the order 
of the Kondo temperature. Substituting Eqn. (|5.3D into 
Eqn. (|5.1|) we obtain the form of the Abrikosov-Suhl 
resonance of the dot system as a function of energy, e: 



P imp (e) 



COs(7T/i) 



1 



where e 



t rkky 2 p cos (2np) + i* + 1 ' (5 ' 4) 

orpRKKY ■ We plot the resonances for two cases, 
2j K 

kdi - e<d > IY2 and \e dl - e d2 \ < r lj2 in Figure 4. 
The plot is for positive energies (as the Bethe ansatz 
solution only directly gives positive energy information) . 
However as we are at the particle-hole symmetric point 
the resonance itself will be symmetric about zero energy. 

In the case \edi — ed2| 3> Ti^, i.e. p ~ 0, we obtain a 
resonance that has the familiar Lorentzian form associ- 
ated with a single doti 36 i 47 (see Figure 4). However this 
resonance has twice the integrated weight in comparison 
with a single dot. Thus even though je^i — e<a| 3> r 12 , 
we cannot think of this resonance as arising solely from 
Kondo physics involving the energetically lower of the 
two dot levels. 

If < |e<ji - ed2\ < ri,2, i.e. n ~ 1/2, the reso- 
nance dramatically changes structure. Instead of peak- 
ing at zero energy it has a sharp peak instead at e = 
2T^ KKY /tt. While the resonance does not vanish at 
zero energy in this case (inset to right panel of Fig- 
ure 4), its value there is dramatically smaller than at 
e = 2T^ KKY /n (by a factor of 1/4cos 2 (tt/!)). The reso- 
nance has thus split about zero energy. This splitting is 
akin to the splitting seen in the Kondo resonance when 
a interdot tunneling term, d\d,2 + d\d\, is turned on be- 

How- 



17.18.26 



tween two degenerate dots, i.e. eai = tdi 
ever unlike this case the splitting is not proportional to 
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FIG. 3: Plots of the T=0 H=0 linear response conductance as a function of gate voltage. In making these plots, we chose 
Ui = 1 and Ti = Ui/20 uniformly and values of \edi — tdi\ to be 5Fi, 3Fi, 2Fi and in panels a, b, c, and d. Y2 and Ui are 
then chosen in accordance with the integrability conditions of Eqn. (|2.8p . On the plots we mark the values of gate voltage 
where the RKKY Kondo effect and the standard Kondo effect take place. 



the bare scale |e<n — e^l, but rather to the much smaller 
dynamical scale, T^ KKY . 

A note of caution is in order. What we have just com- 
puted is the low energy density of states due the impurity 
double dot. But it is not exactly equivalent to the impu- 
rity Green's function 



(4dl)retH + (4 



Formally the two differ by energy dependent matrix el- 
ements (matrix elements which are not directly com- 
putable in the Bcthe ansatz). In the single dot case these 
matrix elements have, at best, only a weak dependence 
on energy. We would guess that this holds in the double 
dot case. But we note that regardless of whether the two 
can be exactly identified (up to a scaling constant), this 
assumption of equivalence does not affect the exactness 
of our T = results for transport . 

We now turn to what sort of Kondo effect is present 
at the particle-hole symmetric point. The Kondo singlet 
that forms at this point involves both electrons sitting on 



the two dots. Because we know that an Abrikosov-Suhl 
resonance forms, we know the electrons are not forming a 
direct singlet with one another (and indeed, in our model 
there is no direct interaction between the dots). Rather 
the singlet formed involving the two localized electrons 
is mediated by the conduction electrons in the leads. We 
thus term this Kondo effect an RKKY Kondo effect. We 
sketch this type of Kondo effect in Figure 5. 

It is clear that the RKKY physics occurring here is 
considerably different than the ferromagnetic interaction 
that fourth order perturbation theory would predict be- 
tween closely spaced dots 4^ We see both that RKKY 
interaction is effectively antifcrromagnetic and occurs at 
an energy scale akin the Kondo temperature. Moreover a 
strong ferromagnetic interaction would be characterized 
by a far different conductance. If the two electrons did 
form an effective spin-1 impurity, the T = conductance 
would instead be 2e 2 /h not the 0e 2 /h we universally find. 

The ground state associated with the RKKY Kondo ef- 
fect is Fermi liquid (as will be obvious when we consider 




FIG. 4: The low lying impurity density of states (the Abrikosov-Suhl resonance) as a function of energy present in the double 
dot system at its particle-hole symmetric point. In the left panel we show the shape (Lorentzian) of the resonance for double 
dots with well separated levels. In the right panel we show the same resonance for the case |e<ji — tdi\ <C ^1,2 (but eji 7^ e^)- 
Here the resonance has split about zero energy by an amount 4T^ KKY /it. 



RKKY Kondo Standard Kondo 




FIG. 5: A cartoon of the two different types of Kondo effects 
that might appear in the double dot: i) an RKKY Kondo 
effect where two electrons total sit on the dots and form a 
singlet mediated through hopping on and off the leads - as this 
effect depends upon edi 7^ td2, the dots possess an unequal 
distribution of electron; and ii) the standard Kondo effect 
where a single electron sits on the two dots localized primarily 
in the lower dot level. 



the low temperature or low field conductances). This 
contrasts with the non-Fermi liquid nature of electrons 
bound into an effective spin-1.— It is a check on our re- 
sults that we do find Fermi liquid physics. We know 
that the properties of a single dot is pcrturbative in the 
Coulomb repulsion and so the dot's low temperature fixed 
point is Fermi liquid. The pcrturbative structure in the 
Coulomb repulsion of a double dot is similar in nature. 
There is then a predilection to believe that the low tem- 
perature physics of the double dot should similarly be 
Fermi liquid. We find further support for our Fermi liq- 
uid picture in Ref. \M- While Ref. studied a two 
impurity Anderson model coupled to two electron chan- 



nels (not one as here), they found that the T=0 electron 
scattering always obeyed the Fricdcl sum rule, a signa- 
ture of Fermi liquid physics, no matter the anisotropy 
between the channels. 

Moving away from the particle-hole symmetric point 
by adjusting edi,2 upwards (keeping their relative dis- 
tance fixed) we uniformly decrease the number of elec- 
trons on the dot until we have in total one localized elec- 
tron. At this point we then expect a standard type Kondo 
effect to take place. This Kondo effect is akin to that 
seen in single level dots at their particle-hole symmetric 
point. While analyzing the Bethe ansatz equations in 
closed form is difficult, we can nonetheless guess what 
the associated Kondo temperature is: 

T Sta„da rd= ^y /2 e ^r^. (5.5) 

This formula, available from the analysis of single level 
dots^&d&SO, is yalid prov ided |e rfl - e d2 \ < Fi,r 2 and 
Ti <^ Ui. In our analysis of the magnctoconductance in 
the next section we verify that indeed ^standard ag gj ven 
above governs the physics. 

When however \edi ~ £d2\ ~ ri,r2, the behaviour of 
the double dot system with one localized electron appears 
Kondo like (as determined by the mangetoconductancc) 
over a much narrower range of magnetic field (in units of 
the Kondo temperature). 

As we further increase the gate voltage, we move into 
a regime characterized by a mixture of interference influ- 
enced by strong correlations on the dots. We see non- 
monotonic behaviour both in the conductance and the 
number of electrons, njis, displaced by the dots. More- 
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over if |e<ji — e^l ^ Ti, T 2 we see that ridis can go nega- 
tive. This requires further explanation. 

The Bethe ansatz does not allow one to direct com- 
pute the number of electrons. rid, occupying the dots. 
Rather it allows one to compute this number. rid, plus 
1/L changes of electron occupancy that arise in coupling 
the leads to the dots. Specifically 



"o-dis = ndol + ndo2 



dx 



c\ a {x)c ea {x)) - p, 



crbulk 



(5.6) 

where p CT buik is the density of electrons of spin a in the 
leads in the absence of any coupling to the dots. In terms 
of transport properties, this is a more natural quantity 
to compute. It is the quantity that is related to the 
scattering phase by the Friedel sum rule, as the proof of 
Langrcth explicitly shows^ We comment on this further 
in the discussion section. 

Now by definition rida, the number of electrons on the 
dots, can never be negative. Thus when we find n CT dis 
to be negative, we know that coupling the dots to the 
leads has induced a negative 1/L correction to the lead 
electron density. 

In this non-monotonic regime we see that sharp varia- 
tions exist in the conductance 2 ^ and are most pronounced 
if kdi — ^d2\ 3> Ti,2- In particular, the variations occur 
at two values of the gate voltage (i.e. tdi)- The first vari- 
ation is associated with a sudden decrease in the number 
of localized electrons n CT( jis &s e ( ii is lowered while the 
second variation sees a sudden increase in n CT( jis- We can 
characterize the two values of edi where these variations 
occur. 

The variations correspond to values of edi where the 
contribution to the scattering phase from Ki 2 is rapidly 
varying (see Eqn. (|4. 1 5[) of Section IV. A). This rapid 
variation occurs when Q, the parameter marking the 
Fermi surface of the A-excitations, is close to JT . Pro- 
vided Q > we can employ Eqn. (|4.16[) to parameterize 
Q in terms of e<ji. Equating Q and l~ x then leads to 
the two values of edi (corresponding to the two values of 
1 1 2) where the variation occurs 



e<n : 



2UiTi 
1 + 



2-rrl: 



1 

2^ 



2tt 



log(2 7 re/r 1 )) 



1/2 



Ui 
' 2 ' 

(5.7) 



The sharpness of these variations in the conductance 
and n^is arc governed by the parameter 

a consequence of K i2 behaving as tan((J~ 1 — Q)/7i) for 
I// 1 — Q small (see Eqn. (|4.15p ). ji vanishes in the limit 
\tdi — 1 — * 00 while becoming large in the opposite 
limit. As the two variations in n,jd\ s arc of the opposite 
sign, the first of the two variations will induce n^is to 
go negative if it does not overlap significantly with the 



second. Thus the condition for n CTC iis to go negative for 
some value of the gate voltage is 

1/r 1 -i 2 ^ \ > 71+72 

that is, the spacing of the variations is larger than the 
sum of their widths. We can translate this into a condi- 
tion upon the necessary separation of the two levels. We 
need 



kdi - £d2 > 



b + \Jb 2 + 4ca 
2a ' 

u?-ui 



b = T 1 T 2 {U 1 + U 2 ) 

c = 2U i Y i \Tl-v\\ 



(5.8) 



if fWis is t° g° negative. If r l5 T 2 <^ Ui,U 2 this condition 
simplifies to 



kdi ~ £d2 > 



8U,T, 



1/2 



(5.9) 



A note of caution is needed here. This condition is de- 
rived under the constraints governing integrability (Eqn. 
(|4.ip ). While self-consistent, it should not be used to de- 
termine the minimum level separation in two dots with 
arbitrary parameters so as to guarantee n^is will be neg- 
ative. However one might nonetheless suppose that this 
condition would provide a rough estimate of the mini- 
mum separation for such dots. 

Up to this point we have discussed in detail the be- 
haviour of the conductance provided edi is not exactly 
equal to td2- If instead e^i = ed 2 , we do not see a marked 
change in the T = 0, H = linear response conductance. 
With edi close to 6^,2, we already sec in Figure 3d) that 
the rapid variations in G associated with interference ef- 
fects disappear. (This is consistent with the slave boson 
study of Ref. (|l9|).) And as edi becomes equal to ed2, 
the functional dependence of G upon the gate voltage 
remains the same (compare panels c) and d) in Figure 
3). However at edi = £d2 there is one significant change. 
At the particle- hole symmetric point, the RKKY Kondo 
effect disappears. As this vanishing takes place discon- 
tinuously, edi = £<i2 marks a first order quantum critical 
point. 

This critical behaviour at edi = £d2 has been ob- 
served before , 33 ' 52 and can be most easily understood by 
first appealing to the non- interacting limit. In the non- 
interacting limit, one can diagonalize the dot degrees of 
freedom via the transformation, 



1 



1/ odd 



(r 1/2 di ±T 2/1 d 2 



and so we have a system where only one dot level is tied to 
the leads' 3 - This mapping is tied to a dramatic change in 
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FIG. 6: The impurity density of states for a non-interacting 
double dot system when tdi and td2 are close but not exactly 
equal and for when €di — tdi ■ We see that at eai — * tdi , there 
is a discontinuous change in the density of states. 

the non-interacting density of states (see Figure 6). We 
see that with e± and ti nearly equal but still different, 
the impurity density of states possesses a sharply peaked 
feature not present when e\ = ti- It is this feature that 
permits the formation of the Abrikosov-Suhl resonance 
when interactions are turned on. Its absence at t\ = e-2 
explains the absence of the RKKY Kondo effect. While 
interactions will couple the second (odd) dot level to the 
system, the absence of direct hopping to the odd level 
means that the RKKY Kondo effect cannot develop. In- 
stead, the two electrons bind directly into a singlet and 
the Kondo effect is entirely absent. 

We point out that the particular behaviour observed 
at £di = e<22 is dependent on there being no coupling 
between the two dots (whether it be tunnel or capacitive 
coupling). If the dots are coupled, the two levels will be 
split by the coupling (for example, a tunneling term will 
lead to split bonding/antibonding dot levels) i 17 i 18 

B. Finite Field Linear Response Conductance at 
the Dots' Particle-Hole Symmetric Point 

We now consider the conductance of the double dots at 
the particle-hole symmetric point. As we have discussed 
in the previous section, this point sees the formation of an 
RKKY Kondo effect with an attendant Abrikosov-Suhl 



resonance. The presence of this resonance governs inti- 
mately how the conductance evolves as Zeeman fields on 
the order of the RKKY Kondo temperature are applied 
to the dots as it represents the sole available low energy 
degrees of freedom. 

The shape of the Abrikosov-Suhl resonance thus deter- 
mines the form of the magnetoconductance as a function 
of magnetic field. In the case that \edi — £d2| 3> IY2, 
the Abrikosov-Suhl resonance has a Lorentzian form (see 
the right panel of Figure 4) with a peak centered at zero 
energy. As a Zeeman field is turned on, the magnetocon- 
ductance will thus immediately depart from its particle- 
hole symmetric value (0e 2 /h). The universal form of the 
G(H) in this case is shown in the rightmost panel of Fig- 
ure 7. We plot both the closed form available from the 
analysis performed in Section IV. B as well as a direct nu- 
merical evaluation of the Bcthe ansatz integral equations. 
We see we obtain an excellent match (thus validating the 
approximations used in solved these equations analyti- 
cally). We see that the magnetoconductance peaks at 
the maximal possible conductance for a value of H cor- 
responding to 1.4Tjf . That is does reach the uni- 
tary maximum is a reflection of the fact that the scatter- 
ing phases for spin up and down electrons are symmetric 
about 2ir (see Eqn. (|4.24[l ). After reaching the unitary 
maximum the magnetoconductance decreases uniformly 
behaving as 1/ log 2 (ff/Tf KKY ) for H > if KKY . 

In the second panel of Figure 7 we plot the magneto- 
conductance for the case je^i — £d2\ <C ^1,2- The func- 
tional form of the magnetoconductance in this case is also 
governed by the shape of the Abrikosov-Suhl resonance. 
As shown in the left panel of Figure 4, the resonance here 
is peaked at an energy e = Tf . While it does not 
vanish at zero energy, its value relative to its peak value 
is 4 cos 2 (tt/j.) where /i is defined in Eqn. (|5.2p and here is 
approximately 1/2. Consequently the change in the con- 
ductance in this case upon the initial introduction of a 
magnetic field will be negligible. But once the magnetic 
field reaches a strength of Tf KKY , the magnetoconduc- 
tance will undergo an extremely rapid variation (as seen 
in Figure 7). Once the unitary maximum is reaches, the 
magnetoconductance again decreases monotonically be- 
having as 1/ log 2 (iJ/f| KKY ). 

We can summarize the behaviour of the magnetocon- 
ductance in these two cases via 



^cos 2 (^) 



H 



o 



H 

jltKKV J , 



if H « T% KKY 



G= < 



21o g 2 (H/T* KKY ) log 4 (,ff/T£ KKY ) 

n ( loglog(H/Tg KKY ) \ (n ( ^loglog(H/rr KY ) \ , n 



(5.10) 



log 5 (ff/T| 



log 6 (_f//Tg KKY 



, ifff»T£ KKY 



r 



where f£ KKY = 2^t£ kky . i n the above expression for the conductance in the limit H ^> Tf KKY we have 
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FIG. 7: Plotted is the magnetoconductance of the double dot system at its particle-hole symmetric point. We consider three 
cases: right panel (\e d i - e d 2\ > Fi^), middle panel (|e d i - e d 2 1 < ri, 2 ), and left panel (|e d i — e d 2\ = 0). 



included the leading order term in /i 2 as this governs the 
different behaviour in the two cases \eai — £<i2\ *C ^1,2 
and \e d i - e d 2\ > I\ 2 . 

We now consider the final case when edi = £d2- As 
discussed previously, in this instance, the Abrikosov-Suhl 
resonance identically vanishes. Without the Abrikosov- 
Suhl resonance and the concomitant absence of a low 
energy density of states, the conductance is unchanged 
upon application of a Zeeman field on the order of T* KKY 
(see rightmost panel of Figure 7). (We note that when 
f-di = £d2 there is no Kondo temperature. The plot in the 
righthand panel of Figure 7 is made taking the Kondo 
temperature as defined in the case |e<ji — e^l ~ but not 
exactly 0.) 



C. Finite Field Linear Response Conductance as a 
Function of Gate Voltage 

We now consider the magnetoconductance away from 
the particle-hole symmetric point. As with this point, the 
changes induced in the conductance by the application of 
a magnetic field is determined by the presence or absence 
of an Abrikosov-Suhl resonance. If a resonance is present, 
we in general marked changes by the introduction of even 
a small Zeeman field. If however the resonance is absent 
(solely in the case |e<ji — e^l = 0), the conductance will 
only be sensitive to Zeeman fields on the order of the 



bare parameters in the model. 

In the top panel of Figure 8 we first consider the case of 
[edi — e<e| S> Ti,2, plotting G as a function of gate voltage 
(i.e. edi) for a variety of values of Zeeman field ranging 
from H = Ti/500 to H — "iTi. As expected from the 
previous section even the smallest of these field strengths 
strongly alters the conductance near the particle hole 
symmetric point (we note that Tf KKY = 0.00243Fi). 
The region of e<ji over which the conductance is altered 
from its zero field value significantly increases as H is in- 
creased. As €di is initially decreased, we sec that the con- 
ductance tracks closely its zero field value (as displayed 
in the top left panel of Figure 3). However a value of edi 
is inevitably reached where the conductance markedly 
changes. For the smallest field value used, H = Ti/500, 
this occurs at edi ~ — 6T1 while for largest, H = STi, 
large scale deviations in G begin at edi ~ Ti- 

For the smaller values of the field used, the conduc- 
tance sees a local maximum around e^i = — 6T1. This is 
the position of the (standard) Kondo effect we argued 
to exist at zero field in the first part of this section. 
The Kondo temperature for this effect is T|- tandard = 
0.00422ri, twice the smallest value of field (H = .002T 1 ) 
displayed in Figure 8. As the magnetic field is increased, 
the conductance in the vicinity of edi decreases as would 
be expected for the standard Kondo effect. 

We plot this decrease in the conductance maximum as 
a function of H in Figure 9. From Rcf. (f36h . the expected 
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FIG. 8: Plots of the conductance as a function of tdi for a variety of values of Zeeman field, H. We consider three cases: a) 
\edi — fd2| > Ti.2; b) \eai — e<j2| -C Ti,2; c) \cdi — ed2 ] = 0. In all three Ui = 1 and Ti = [7/20. For a) and b) we see sharp 
changes in G closely tied to the presence of an Abrikosov-Suhl resonance at H — 0. For c) the resonance is absent and changes 
in G occur only when H exceeds Ti. 
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FIG. 9: A plot of the magnetoconductance for the case \edi — 
tdi\ 3> ^1,2 at a value of the gate voltage where one localized 
electron is found on the dots. We use the same parameters as 
in Figure 8. We plot both the results of a numerical analysis of 
the Bethe ansatz equations as well as the small field analytic 
form that would be expected if the system was a single dot. 

from of this decrease is 

G(H) = G(H = 0)(1 - ^ (jslasi) \ ( 5 ' n ) 

We see in Figure 9 that the numerical analysis of the 
Bethe ansatz equations follows this form for values of 
the Zeeman field up to Tf tandard . 

The plot in Figure 9 was produced in the case je^i — 
e d2\ ^ r 12 . As this condition is relaxed, the standard 
Kondo effect does not disappear. However the region of 
field over which Eqn. (|5.1ip is valid decreases. 

In the middle panel of Figure 8 we consider the magne- 
toconductance in the case \edi — tdi\ <C In this case 
the conductance behaves in much the same fashion as 
\^d\ — 1 ^ Ti,2- For the smallest values of the field the 
conductance goes unaltered from its zero field value over 
a substantial range of gate voltages. This range shrinks as 
H increases. And when the conductance ceases to track 
its H = value, it does so abruptly. Due to the nature of 
the Abrikosov-Suhl in this case (i.e. a resonance sharply 
peaked about e = T^ KKY ), we find that the changes in 
the linear response conductance, when they occur, to be 
even more marked. Even for fields for in excess 

ofT RKKY 

(in this case T^ KKY = 0.001 79Fi), we see sharp changes 
persist in the magnetoconductance. For both H = Y\ 
and H = ?>Ti we see that there are sharp downward 
changes in the conductance for values of e d i near 0. 

In the final panel of Figure 8 we plot the magnetocon- 
ductance in the case \e d i — e d 2\ = 0. Here there is no 
Abrikosov-Suhl resonance at the particle-hole symmetric 
point. With the absence of a dynamically enhanced low 
energy density of states, we find G to be insensitive to 
values of H smaller than the bare parameters. The con- 
ductance is approximately equal to its zero field value 



until H becomes of the same magnitude as T±. For val- 
ues of H in excess of Ti we see that the zero field peak 
at tdi ~ — Ti) splits in two with increasing separation as 
H grows. 

VI. FINITE TEMPERATURE TRANSPORT OF 
DOUBLE DOT 

In this section we analyze the finite temperature linear 
response conductance, G(T), at the particle- hole sym- 
metric point. As with the magnetoconductance at this 
point, the behaviour of G(T) at temperatures on the or- 
der of T^ KKY is ultimately determined by the form of 
the low energy density of states, that is, the Abrikosov- 
Suhl resonance as given in Eqn. (|5.4p . In particular, the 
nature of the Abrikosov-Suhl resonance marks the differ- 
ence between the two opposing limits of the bare level 
separation, i.e. \e d i - e d2 \ > T 1:2 and |e dl - e d2 \ < T ia . 

To evaluate the finite temperature conductance, we 
solve numerically through iteration the equations in 
(|3.22p and (|3.26|) . To do so we approximate the source 
term for the integral equation governing p lmp (fc) by Eqn. 
(|5.1[) , an excellent approximation provided we are at the 
particle-hole symmetric point. Then following the pro- 
cedure described in Section III.C, we obtain the results 
plotted in Figure 9. 

In the case \e d i — £ d2 \ ^ (see the left panel of 
Figure 9), the Abrikosov-Suhl resonance is peaked at zero 
energy Thus even for temperatures far below the Kondo 
temperature, T^ KKY , we see an immediate change in the 
conductance from it zero temperature value (0e 2 /h). For 
temperatures far below K the conductance has a 
quadratic Fermi liquid form which we can determine to 
be 

= — 4 (,2*887 ) ■ (6-1) 

We can determine this form by comparing the analy- 
sis of the Bethe ansatz equations at finite temperature 
for a single level dot.— For a single level dot (sld) the 
Abrikosov-Suhl resonance is given by 

p i m psld (e! ^d )= 1 1__ (6>2) 



For a double dot (dd) with well-separated bare levels (i.e. 
/1 ~ 0), the Abrikosov-Suhl resonance has the same form 
(see Eqn. (pT4")) ) 

p imp dd (£ ^RKKY) = 2cos( M 7T)p im P sld (e, T« KKY ). 

(6.3) 

This in turn implies that an electron with energy e scat- 
ters off the double dot via the relation 

<5 d 1 d (T, e ) = 2cos( / i7r)^ d (T, e ). (6.4) 

This relation follows because the Abrikosov-Suhl reso- 
nance serves as the source term governing the linear in- 
tegral equations determining the scattering phase. For a 
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FIG. 10: Plots of the finite temperature linear response conductance in the cases \e<ii — ta,i\ ^> ^1,2 (left panel (a)) and 
|e<ji — £d2| <C Ti,2- In both cases we employ Ui = 1 and Fi = U/20. For (a) we use a bare level separation of ed2 — £di = 5Fi 
while for (b) ed2 — tax = 0.5Fi. 



single level dot, the deviation of the finite temperature 
linear response conductance from its zero temperature 
value takes the fori m 36 ' 50 



5G(T) 



16 (T s K ld ) 2 ' 



As G(T) = J d E f(E) sin 2 (6 c] (E,T)), we see that for 
6 \(E) small, the low temperature form of G(T) for the 
single level dot determines G(T) for the double dot to be 
given as in Eqn. (|6.ip . 

As the temperature is increased, G(T) reaches a max- 
imum at roughly T = T^ KKY . With further increases, 
G(T) decreases from its maximal value (for the parame- 
ters chosen, lAe 2 /h) gradually. For large values of tem- 
perature we see that the conductance decreases logarith- 
mically. From the behaviour of G(T> Tk) for a single 
level dot (Ref. (|54l)) and the similarity of analysis be- 
tween that case and a double dot, we might expect 



G{T > T* 



RKKYn 



RKKY ^ 



(6.5) 



However from examining panel (a) of Figure 9, it is 
clear that if this behaviour docs hold, it is valid only for 
temperatures much larger than 10T^ KKY (in the range 
10T rkky_ 100r RKKY ; Figure 9 i mr> li e s the conductance 

is decreasing at best according to 1/ log(T/T^ KKY )). 
This however is unsurprising as a similar condition holds 
on the finite temperature behaviour of a single level dot 
for T~>T K . 

In panel (b) of Figure 9 we plot G(T) for the case 
of |edi — €^2 1 <C T\,2- In this case the Abrikosov-Suhl 
resonance has a sharp peak at an energy e = T^ KKY with 



minimal (although finite) spectral weight at zero energy. 
We thus expect that the low temperature response to be 
much smaller than for the case |edi — e<22| 3> Ti.2- And we 
see this in Figure 9 where G(T) is effectively Oe 2 /h (equal 
to its zero temperature value) until T ~ 0.1T^ KKY . For 
T < T^ KKY , G(T) has the exact same form as given 
in Eqn. (|6.1[) for the case |e<ji — ed2| 3> ^1,2- However 
here n ~ 1/2 and so the coefficient of the T 2 correction 
is extremely small. 

As the temperature is increased in this case to roughly 
2T^ KKY , a maximum occurs in the conductance (for the 
particular parameters chosen). While the value at which 
the maximum occurs is non-universal, we do expect it 
to roughly occur when the temperature is 0(T^ KKY ). 
With further increases in temperature, the conductance 
logarithmically decreases. If this decrease is to behave 
as G ~ l/log 2 (T/T^ KKY ) (as for a single level dot), the 
temperature at which this asymptotic form begins to be 
valid is far in excess of 10T^ KKY . 



VII. DISCUSSION 

A. Friedel Sum Rule 

The Friedel sum rule relates the scattering phase, 5 ae , 
of an electron of spin a at the Fermi surface at zero tem- 
perature to the number of electrons displaced by connect- 
ing the leads to the dot, n a ^ss\ 



S<Tt 



2mi n 



(7.1) 
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n CT dis has two contributions: 1) the number of localized 
electrons sitting on the dot and equal to, 



and 2) the l/L deviations induced in the electron number 
in the leads 



Pleads = dx 



where Po-buik = L/2it is the bulk interacting density of 
states of the leads with the leads and dots unconnected to 
one another. While this second contribution is typically 
not considered in the application of the Friedel sum rule, 
its proof explicitly takes it into account.— 

We have already seen that in the model of double dots 
analyzed here that this second contribution, <5ni ca d s , is 
non-negligible and in fact can be negative, driving TVdis 
itself negative (and thus cannot be ignored if the conduc- 
tance is to be computed correctly). The range of gate 
voltages where this occurs marks a region of rapid varia- 
tion in the conductance (see the top left panel of Figure 
3 for example). The possibility that this second contri- 
bution is non-zero is closely tied to the presence of two 
interfering pathways for the electrons to transit the dot 
system. A finite value for Sn\ ca ds has also been seen in the 
analysis of Fano resonances in quantum dots.— For dots 
exhibiting Fano resonances, there are two tunneling paths 
between the leads: one resonant inducing a single parti- 
cle scattering phase, 5(q) = — 2tan~ 1 (r/(g — e^)) and 
a non-resonant second path with a constant scattering 
phase. The two paths, however, interfere and the overall 
scattering phase is more complicated than the mere sum 
of the two (akin to the fact the sum of the bare scattering 
phase off the double dots in Eqn. (|2.4j) has the sum in 
the argument of tan -1 ). The system exhibiting Fano res- 
onances may be considered a certain limit of the double 
dots studied here. If we take the limit of tdi and T2 in 
such a way that 



To 



q - £d2 



is constant for all wavevectors q we recover the physics 
of Fano resonances. 

TVdis is a quantity naturally available from the Bcthe 
ansatz as the Bethe ansatz allows one to directly access 
the scattering phase, i.e. the formalism, as discussed in 
Section III, provides a means to derive the scattering 
phase independent of directly invoking the Friedel sum 
rule. In this sense we provide a proof of the Friedel sum 
rule in the case of a double dot system. We do note that 
however we are able to compute in principle the occupa- 
tion of the dots YlaWradaa) by computing the impurity 
free energy and then taking appropriate derivatives with 
respect to e ( i a . 



B. The Effects of Breaking Integrability 

We now turn to the question of how fined tuned the 
results presented herein are. Our ability to exactly solve 
the models of the double dots depends upon the con- 
straints listed in Eqn. (|2.8[) being satisfied. However 
it is important to determine how the physics is affected 
if these conditions are broken, at least weakly. While 
it is easily conceivably that an experimental double dot 
system could be tuned so that its parameters satisfy ap- 
proximately the constraints of Eqn. (|2.8[) . we cannot 
hope that the constraints are met exactly. 

Fortuitously, we do not expect any of the major phe- 
nomena identified in this paper to change qualitatively 
if Eqn. (|2.8p is weakly violated. Generally, we expect 
the physics to be perturbative in any deviation from the 
conditions necessary for integrability. To see this, con- 
trast the situation of a single dot coupled to a single 
lead. This problem is non-perturbative in the dot-lead 
coupling, a partial consequence of the ambiguous nature 
of the zeroth order zero temperature dot density ma- 
trix about which one conducts the perturbation theory. 
This is true both in the presence and absence of inter- 
actions. Indeed it is relatively straightforward to derive 
the dots' Green functions at U = and see that they 
do not have a well-defined V — > limit. Concomitantly, 
perturbation theory in U , the dot Coulomb repulsion, is 
well behaved as the dot Green's function is unambigu- 
ously defined at U = by the presence of the dot-lead 
hybridization. Similarly if we weakly break one of the 
integrability constraints, we expect the deviation to be- 
have in a perturbative fashion as the dot state is already 
robustly established. This principle guided the study of 
perturbations about the exactly solvable Toulouse limit 
of a non-equilibrium dot in Rcf. HH). 

Now while small deviations from the integrability con- 
straints will only quantitatively change the physics, the 
same cannot be said of large deviations. In Ref. (|28f ) a 
particle-hole symmetric double dot is considered where 
(-di = —Ui/2, 6d2 = U2 = 0, and Vi, V 2 finite. These con- 
ditions represent a strong violation of the integrability 
condition, V 2 2 U2 = V 2 U\- And instead of a Kondo effect, 
the double dot is found to be in a local moment regime. 

There are other more specific arguments as to why vi- 
olating the integrability constraints will not drastically 
change the physics. To this end, let us consider the con- 
straint upon the left-right hopping amplitudes. If this 
constraint is met, only one electron channel couples to 
the dot (the even channel), while the second channel is 
decoupled. Let us consider the affects of weakly violating 
this condition. The second (odd) channel then couples to 
the dots. However this will change qualitatively none of 
the Kondo physics discussed in the manuscript. In prin- 
ciple, the second channel could introduce a second Kondo 
scale, Tk2 (for either the RKKY Kondo effect or the reg- 
ular Kondo effect) as in Rcf. (f56h . If we start at some 
high temperature and begin to lower it, when T reaches 
Tki (the Kondo temperature associated with the even 
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channel and given in the manuscript) the Kondo singlets 
(either RKKY or standard) will establish themselves. If 
we continue to lower the temperature to Tki nothing fur- 
ther will happen as there are no spin degrees of freedom 
(the state is a singlet) left to screen. (This contrasts with 
Ref. jHB) where an effective spin-S (S > 1) impurity was 
envisioned and a two-stage Kondo effect was predicted to 
occur). 

Another argument that demonstrates the RKKY 
Kondo effect is not a result of fine tuning makes ap- 
peal to a Schrieffcr- Wolfe transformation at the particle- 
hole symmetric point. Under such a transformation the 
Hamiltonian reduces to the form, 

H = -ffcond uctionclcctrons 

where the two couplings, Ji, are given by 



Ji 



2r, _ 2T* 
U + e di e di 



We see that generically J\ ^ J2 for the intcgrablc range 
of parameters and that deviations away from the inte- 
grability constraints do not change this fact. Thus we 
would expect such deviations to affect only non-universal 
aspects of the physics (such as the exact values of Kondo 
temperatures). 

While a small breaking of the integrability conditions 
will not destroy the RKKY Kondo effect, there is one 
small perturbation that will. The RKKY Kondo effect 
involves singlet formation of two electrons sitting on two 
dots as mediated by the electrons sitting in the leads. 
But if one were to add a coupling, J s i n giot, in excess of 
the Kondo temperature, T^ KKY , that promoted direct 
singlet formation between the two localized electrons, i.e 
SH = Jsingiot £1 ■ S%, we would expect to destroy the 
RKKY Kondo effect. In principle, any experimental real- 
ization of a double dot system will contain such a J s i n giet- 
And because T^ KKY is exponentially small compared to 
the bare parameters, one might think that the RKKY 
Kondo effect will be generically overcome. While this 
may be true at the particle- hole symmetric point, the 
RKKY Kondo effect exists over a range of gate voltages 
centered around this point. By tuning the gate voltage 
to different values within this range, one can increase 
T« KKY above any putative value of Jsingiet- 

We finally turn to the effects of breaking the integra- 
bility conditions on the the quantum critical point at 
(-di = £d2- We note that this point is tied to the dis- 
appearance of a singular feature in the non-interacting 
impurity density of states at edi = ^di (sec Figure 6). 
In turn this disappearance is dependent upon the de- 
coupling of the odd dot degree of freedom, d a dd- A 
perturbation that restores this coupling (such as weakly 
coupling a second channel of electrons to the dots) will 
change the critical point to a smooth crossover and so re- 
store the RKKY Kondo effect, but a perturbation which 
leaves d dd uncoupled (such as perturbations on the inter- 
/intra-dot Coulomb repulsion) will presumably leave the 
critical point in place. 
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APPENDIX A: WIENER-HOPF ANALYSIS OF 
EQUATIONS GOVERNING THE LINEAR 
RESPONSE CONDUCTANCE 

1. Review of Wiener-Hopf Analysis 

We review the method of solution as presented in Ref. 
35h for equations of the general form 



/(*) 



dz'f{z')h(z-z')+g(z). (Al) 



Writing f ± (z) = f(z)8(±z =F A), the Fourier transform 
of the above equation yields 

+ r M = f-[u)h(u) + g(oj), (A2) 

where Fourier transforms are defined by 



0(0;) = / du;e iuJz a(z). 



(A3) 



The key step in the analysis is writing 1 — h(uj) as a 
product of functions, G ± , which are analytic in the up- 
per/lower planes respectively: 



1 - h(w) 



1 



G+(w)G-(w) 
We can then write Eqn. (|A2[ ) as 



(A4) 



G (LO) 

(A5) 

Given e lu ' A f ± (tjj) is analytic in the upper/lower half 
plane, applying the operators 



2ni J u>' — (to ± iS) ' 
to Eqn. (|A5|) yields solutions for both /+ and /" 



(A6) 



/"(to) = 



x / dw' 



1 



2?ri 



lj' — (lu — iS) 
e iuA 1 



5 ( W ')G + (^')e 



/\„— iui'A. 



x / dtu' 



' i\r<+t. ,'\„—iu'A 



oj' — (lu + iS) 



ff ( W ')G + ( W ')e^ A . (A7) 
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2. A Wiener-Hopf Analysis of the Scattering Phase 
Analysis of Zero Field 

The scattering phase can be written as 
6 e = 2tt / dA(<ri(A) + (72(A)); 



= 27r-7r / d\(ai(X) +<r 2 (A)). (A8) 

J — oo 

where each <xs(A) satisfies 



a™ p (A) = / dX'R(X- A> ? imp (A') 

— OC 

/oo 
dA'(l + a 2 )- 1 (A-A')a |1 _ /3 - 1| (A') 
-OO 

+ / dgs(A-5(9))Aoi(g). (A9) 



The above two equations can be obtained by inverting 
via Fourier transform Eqn. (|4.8p . 
Identifying 



i 



1 + e-H 



e" M/2 / d g A 0l (g)e M9( «) 



+sign(l-/3i)- 
1 



1 + e~H 



-(l-Z^signa-A); 



2t/ i r i 



2*7T 



r(l/2=F iuj/2n) \ 



, (A10) 



where Pi and cti are defined in Eqns. (|4.5p . wc apply the 
Wiener-Hopf methodology presented at the beginning of 
this section to compute 



dXai(X) = o-r( w = 0). (All) 



Wc so find 



where 



rQ 

/ dXai(X)=K a +K i2l (A12) 

J — OO 



xr(- + ^) / d q & 0l { q y u{9{q] - Q) 

2 27T J_ 00 

1 /• 1 

K i2 = sign(/3, — I) , / — — f (-£cJ + 5) 



2w 3 / 2 i J lo + id 
-2 27r n 27re' 



x(iw + *) - r(- + —)(—) **e luj{Q '« >. 



(A13) 

We will evaluate each term K a / 2 in turn. 

For A'ii, if Q < 0, we can continue the contour into 
the upper half plane with the result 



K,(Q <0 ) If til! (ZLtVW/a 

/OO 
dgA 0l ( g )e-"( 2n+1 ) (9 ^- Q ). (A14) 
-oo 

If instead Q > 0, we directly evaluate the integral J dq, 
dqA 0l (q)e luJ9 ^ « e^ 1 * -1 "^ M . (A15) 



Setting Ji = I i 1 — Q we can then write if,i as 

K a = -^r- / — w — r - nr . (A16) 



V2tt 



G- (w)(w + te) 



If Ji < 0, we evaluate this integral by continuing into the 
lower half plane: 

1 f°° 1 
K a (Q > 0, J, < 0) = 1 - J duT(- + u) 

^ sin(7T^(l + P' 1 ) e y &2vJiU) 



= l -^§r +oi ^ < A17) 

where in the last line we have written down an expansion 
in Ji to leading order in 1/J*. If instead J; > 0, i.e. 
< Q < I^i we close the contour in the upper-half 
plane obtaining, 



dio sin(7r/3,- lj) fu)\ e 



i io w cos(7ro;) \ejT(^+uj) 
1 ~ (_i)» /- ■ 

. ^-^ nil 



^ cos(^ri(n+l)) e " 2 ^("+^) 



r(l + n) 

where the symbol P indicates that the principal value of 
the integral is to be taken. We can alternatively evaluate 
this integral and sum in a fashion similar to that used in 
Refs. (HI) and HI) to obtain, 

K ll (Q>Q,J l >0) = y/2-1- -^Jt + i^J? 
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f 



+- 7f - : (/3r 1 - m + om 1 - i) 2 Jt) + o{jf). (A19) 

We similarly evaluate if ,2 for Ji > and J,; < 0. For 
Ji > wc can continue into the upper half plane with 
the result 

K a (Ji>0) - Sign(1 ' A) P Hm-)"***" 
V* Jo e 
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sin(27rw(7i - |)) 



T(i + w) cos(7toj) 
sign(l - A) J2 



n— 2 



1 \ n+» 
2 



- 27rJ ' w cos(7r(2n + - §) 



r(l + n) 
sign(l-fl)( 7i -l/2) 



+ (A20) 



'(g) / dA'i?( ff ( 9 ')-ff(9))pr P (9')- (A24) 



As we are interested in the behaviour of /5™ p ((?) at q < 0, 
we introduce the change of variables z = g(q) for 2 > 



and q < and define p" np (z) ; 
above equation is then recast as 



-P l n<l)l9'{<l)- The 



And if Ji < 0, we instead continue into the lower half 
plane with the result, 



K i2 (Ji<0) = sign(l-A) 1 



1 



2-nJiUj 



duoe 



LO 2 £J 

= sign(l-A)(l-^)+0(ij)). (A21) 

We now finally turn to how Q depends on f/j, and 
e<jj. As the bulk density equations are the same as for 
a single dot, this dependence can be taken directly from 
Refs. |35l ) and If Q < 0, Q is determined implicitly 
by the equation 

2^ 777T~TTvT79 G (iir(2n + 1)). 



vWI 7T ^ (2n + 1)3/ 2 



If further Q ^> 0, Q is given by 

1 



(A22) 



Q = q* + 77- ln(27req* 

Z7T 



(A23) 



We note that there are no inconsistencies in the above 
equations as the the quantities 2UiTi and Ui + 2t d i are 
independent of i = 1, 2. 



3. Analysis of Finite Field Equations 

The equation we must solve takes the form 

imp/ \ imp/ \ 

Pi W = Poi W 



imp / \ 

Pt ( z ) 



imp / \ 

p 0i ( z > 



dzR(z- z')p imp {z'), (A25) 



where b = B 2 /(2UiTi) and we can write the Fourier 
transform of p^ p (z) as 



imp / \ 

Pm ( w ) 



due^p^{z) 

cosh( W (| - ^-))e^ 



e(A - 1)- 



cosh(f) 



-|tanh| / 



The kernel, R{z) 1 of the integral equation in Eqn. (|A25j) 
is the same encountered in analyzing the scattering phase 
at zero field. Thus using the analysis at the beginning 
of this Appendix together with the properties of R{z) 
already established, we can immediately write down the 
form of p 1 ^ ap (z): 



dqpT P {<l) - 



G+(0) 
2?ri 



du> 



id 



(A27) 



Substituting in the above form for p™ p (a;), we obtain 
the solution as written in Eqn. ()4.28|) . As discussed 
in Section V.B, the Fermi surface parameter, &, is the 
same as found in the analysis of the Anderson model 
for single level dots and for small fields takes the form 
b = — - log 7reg2 
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